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1.0  SCOPE 

1.1  Introduction 

The  Global  Ocean  Forecast  System  (GOFS)  Version  3.0  (V3.0)  consists  of  the  1/12° 
global  HYbrid  Coordinate  Ocean  Model  (FIYCOM)  and  uses  the  Navy  Coupled  Ocean 
Data  Assimilation  (NCODA).  This  User’s  Manual  (UM)  describes  the  sequence  of  scripts 
run  daily  in  the  normal  processing  stream  at  the  Naval  Oceanographic  Office 
(NAVOCEANO)  and  is  geared  toward  those  users.  A  more  generalized  HYCOM  User’s 
Manual  for  running  simulations  in  research  mode  may  be  found  on  the  HYCOM 
consortium  website  (http://hycom.rsmas.miami.edu).  Accompanying  documents  to  this 
manual  are  the  Software  Design  Document  (SDD)  (Wallcraft  et  ah,  2009)  and  the  Phase  I 
Validation  Test  Report  (VTR)  (Metzger  et  ah,  2008).  It  should  be  noted  that  when  the 
VTR  was  written,  the  system  as  a  whole  was  referred  to  as  Global  Ocean  Prediction 
System  (GOPS)  V3.0.  However,  NRL-developed  global  ocean  nowcast/forecast  systems 
have  since  been  renamed  the  Global  Ocean  Forecast  System  (GOFS).  It  is  referred  to  as 
such  in  the  SDD  and  this  UM.  However,  GOPS  and  GOFS  are  interchangeable  between 
these  documents. 

The  GOFS  V3.0  runstream  is  depicted  in  Figure  1.  Relative  to  the  nowcast  time,  the  first 
NCODA  ocean  analysis  is  perfonned  at  t  =  -126  hours  with  the  analysis  window  for 
altimeter  (and  all  other)  data  spanning  ±36  (±12)  hours.  (The  first  hindcast  goes  back  5+ 
days  from  the  nowcast  due  to  late-arriving  satellite  altimeter  data.  An  examination  of  the 
timeliness  of  the  historical  altimeter  data  determined  an  additional  data  gain  of  18% 
between  four  and  five  days;  orbits  also  improve  with  the  age  of  the  data.). 


Manuscript  approved  February  4,  2009. 
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NCODA  analysis  windows  centered  at  these  times 
±36  hours  for  altimeter  data 
±1 2  hours  for  all  other  data 

1)  Perform  first  NCODA  analysis  centered  on  tau  =  -126 

2)  Run  HYCOM  for  24  hours  using  incremental  updating  (□ )  over  the  first  6  hrs 

3)  Repeat  steps  1 )  and  2)  until  the  nowcast  time 

4)  Run  HYCOM  in  forecast  mode  out  to  tau  =  120 


Figure  1:  The  HYCOM/NCODA  runstream.  Approximate  run  times  using  379  IBM 
Power  5+  processors:  a)  six  NCODA  analyses  -  1.1  hours/analysis  =  6.6  hours,  b)  five 
HYCOM  hindcast  days  using  a  240  second  timestep  -  0.8  hours/model  day  =  4.0 
hours,  c)  five  HYCOM  forecast  days  using  a  240  second  timestep  -  0.8  hours/model 
day  =  4.0  hours,  for  a  total  of  d)  14.6  wall-time  hours. 


After  the  initial  NCODA  analysis,  HYCOM  is  run  for  24  model  hours,  with  the  NCODA 
analysis  incrementally  updating  the  ocean  model  over  the  first  six  hours.  Therefore,  at 
00Z,  HYCOM  has  fully  ingested  all  observational  data  used  in  the  analysis.  The  NCODA 
analysis  and  HYCOM  hindcast  cycle  repeat  themselves  daily  up  to  the  nowcast  time. 
HYCOM  then  continues  to  run  in  forecast  mode  out  to  120  hours  (5  days).  The  forecast 
length  of  five  days  was  chosen  in  part  because  atmospheric  forcing  typically  is  available 
that  far  out  in  time.  As  currently  configured,  the  atmospheric  forcing  is  blended  toward 
climatology  after  the  forecast  output  ends,  so  the  system  could  run  in  forecast  mode  out  to 
-30  days.  However,  this  would  require  considerable  additional  computational  resources. 


2.0  REFERENCED  DOCUMENTS 

2.1  Government  and  Technical  References 

Cummings,  J.  and  Carroll,  S.,  (2006).  Software  User's  Manual  for  the  Navy  Coupled 
Ocean  Data  Assimilation  (NCODA)  System,  NRL  Tech.  Rpt.,  MRY-001-06,  Naval 
Research  Laboratory,  Monterey,  CA. 

Metzger,  E.J.,  Smedstad,  O.M.,  Thoppil,  P.,  Hurlburt,  H.E.,  Wallcraft,  A.J.,  Franklin,  D.S., 
Shriver,  J.F.  and  Smedstad,  L.F.,  (2008).  Validation  Test  Report  for  the  Global 
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Ocean  Prediction  System  V3.0  -  1/12°  HYCOM/NCODA:  Phase  I,  NRL  Memo. 
Rpt.,  NRL/MR/7320 — 08-9148,  Ocean  Modeling  Division,  Naval  Research 
Laboratory,  Stennis  Space  Center,  MS. 

Wallcraft,  A.J.,  Metzger  E.J.  and  Carroll  S.N.,  (2009).  Software  Design  Description  for  the 
HYbrid  Coordinate  Ocean  Model  (HYCOM)  Version  2.2.  NRL  Memo.  Rpt., 
NRL/MR/7320 — 09-9166,  Ocean  Modeling  Division,  Naval  Research  Laboratory, 
Stennis  Space  Center,  MS. 
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3.0  APPLICATION 

3.1  Description  of  GOFS  Usage 

This  manual  describes  the  procedures  for  running  the  Global  Ocean  Forecast  System 
Version  3.0.  It  is  designed  for  NAVOCEANO  personnel  who  will  be  running  GOFS  in 
their  operational  queues.  The  user  may  also  refer  as  needed  to  the  HYCOM  Software 
Design  Description  (SDD),  which  describes  the  code,  physics  and  basic  equations 
(Wallcraft  et  ah,  2009),  and  the  Validation  Test  Report  (Metzger  et  al.,  2008)  for  this 
system.  Additionally,  the  NCODA  UM  (Cummings  and  Carroll,  2006)  may  be  of  use. 

3.2  Memory  and  Processor  Allocation 

HYCOM  can  run  on  any  number  of  processors  and  its  footprint  can  be  changed 
accordingly.  It  is  currently  configured  to  use  619  processors,  but  596  or  713  processors 
are  other  multiples  for  which  the  code  is  compiled.  Memory  requirements  for  NCODA 
vary  depending  upon  the  amount  of  data  that  goes  into  the  analysis. 

3.3  GOFS  V3.0  Runstream  Scripts 

For  a  daily  integration  of  GOFS  V3.0,  a  sequence  of  scripts  are  generated  and  run.  All 
scripts  contain  environmental  variables  ${E}  and  ${EXPT}  which  define  the  experiment 
number.  The  example  that  follows  uses 

setenv  E  727 

setenv  EXPT  expt_72.7. 

If  a  new  experiment  number  is  desired,  these  variables  would  have  to  be  changed  in  all 
scripts  and  the  appropriate  directory  structures  would  need  to  be  created.  Figure  2  is  a 
schematic  flowchart  of  the  system  and  all  scripts  are  provided  in  Appendices  A-F. 

Several  scripts  are  used  to  direct  the  running  of  GOFS  V3.0.  These  are  summarized  in 
Table  1  below.  The  script  names  often  begin  with  the  experiment  number,  in  this  example, 
727. 
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Table  1:  GOFS  V3.0  .csh  and  .com  scripts. 


Script _ _ Description 


727shellsub_submit_ncoda.csh 

Main  driver  script.  Generates  scripts  for  each  daily  run.  See 
Section  3.3.1  and  Appendix  A. 

7271sf_wind_prep_$  {idtgtod}  1 8 
$  {idtg}  18.com 

Atmospheric  forcing  generation  script.  See  Section  3.3.2 
and  Appendix  B. 

727extr_${subreg}_$  {idtgtod}  18 
_$  {idtg}  18. csh 

Extracts  NCODA  z-levels  and  surface  fields  from  the 
HYCOM  archive  file  for  each  NCODA  subregion.  See 
Section  3.3.3  and  Appendix  C. 

727ncoda_$  {subreg}_$  {idtgtod} 
18  $  {idtg}  18. csh 

Runs  the  NCODA  analysis  for  each  subregion. 

See  Section  3.3.4  and  Appendix  D. 

727ncoda2arch_$  {idtgtod}  18_${ 
idtg}  18.  csh 

Converts  the  NCODA  analysis  on  z-levels  to  HYCOM 
vertical  coordinates.  See  Section  3.3.5  and  Appendix  E. 

7271sf_ncoda_$  { idtgtod}  18_${  id 
tg}  18.com 

Runs  HYCOM.  See  Section  3.3.6  and  Appendix  F. 

3.3.1  Main  Driver  Script 

The  main  driver  script  may  be  found  at 

/u/home/ooc/models/hycom/GLBa0.08/expt_72.7/mvoi/727shellsub_submit_ncoda.csh. 
Running  as  a  cron  job,  it  generates  the  set  of  scripts  needed  for  each  daily  run.  Most  of  the 
work  is  simply  script  creation,  but  its  final  task  is  submitting  a  job  to  create  the 
atmospheric  forcing  needed  by  GOFS  V3.0.  A  subsequent  sequence  of  scripts  then  gets 
executed.  Appendix  A  provides  the  complete  script  for  viewing. 


The  following  command  sets  the  number  of  hindcast  days.  It  is  currently  set  to  five  days. 


# 

#  Go  back  5  days  from  today  to  do  a  5  day  hindcast 

#  idtg:  date-time-group  for  today  minus  5  days 

# 

set  idtg= ' /u/home/$ { user } /bin/addndays  yyyymmdd  ${idtgtod}  -5' 
# 


3.3.2  Atmospheric  Forcing  Generation  Script 

The  atmospheric  forcing  generation  script  is  found  at 

/scr/ooc/data/hycom/GLBa0.08/expt_72.7/logs/7271sf_wind_prep_${idtgtod}  1 8_$  {idtg}  1 
8.  com 


5 


NRL/MR/7320— 09-9175 


GOFS  V3.0  User’s  Manual 


Here  environmental  variable  $  {idtgtod}  is  the  date-time-group  for  today  (e.g.  20081206) 
and  variable  ${idtg}  is  the  date-time-group  for  today  minus  the  hindcast  length  (e.g. 
20081201).  See  Appendix  B  for  a  copy  of  the  entire  script. 

3.3.3  NCODA  Z-level  Extraction  Script 

This  script  extracts  NCODA  z-levels  from  the  HYCOM  archive  file. 
/scr/ooc/data/hycom/GLBa0.08/expt_72.7/logs/72.7extr_${subreg}_$  {idtgtod}  18_${idtg 
}18.csh 

Environmental  variable  ${subreg}  is  the  subregion  identifier.  Currently  there  are  seven: 
ANTarc,  ARCatl,  ARCocn,  ARCpac,  MERatl,  MERinl,  and  MERpal.  Appendix  C 
provides  the  complete  script  and  an  example  from  MERpal. 

3.3.4  NCODA  Subregion  Run  Script 

This  script  runs  the  NCODA  analysis  for  each  subregion. 
/scr/ooc/data/hycom/GLBa0.08/expt_72.7/logs/727ncoda_${subreg}_$  {idtgtod}  1 8_$  {idtg 
}  18.csh.  All  subregions  are  run  in  parallel.  These  jobs  also  start  the  NCODA  plotting  and 
file  transfer  to  the  permanent  storage  machine.  See  Appendix  D  for  the  complete  script  and 
an  example  from  MERpal. 

3.3.5  NCODA  Z-level  to  HYCOM  Vertical  Coordinate  Conversion  Script 

This  script  converts  the  NCODA  analysis  on  z-levels  to  HYCOM  vertical  coordinates. 
/scr/ooc/data/hycom/GLBa0.08/expt_72.7/logs/727ncoda2arch_$  {idtgtod}  18_$  {idtg}  1 8. c 
sh.  This  job  is  submitted  by  one  of  the  NCODA  analysis  scripts  and  will  wait  until  all  the 
subregions  are  completed  before  it  executes.  Upon  completion,  it  submits  a  plotting  job  to 
HYCOM.  Appendix  E  shows  a  copy  of  this  script. 

3.3.6  HYCOM  Run  Script 

/scr/ooc/data/hycom/GLBa0.08/expt_72.7/logs/7271sf_ncoda_$  {idtgtod}  18_$  {idtg}  18. co 
m  is  the  script  that  runs  HYCOM.  At  the  end  of  this  script,  the  next  day  is  submitted  by 
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727shellsub_daily.csh,  which  is  very  similar  to  the  cron  script  that  initiates  the  entire 
sequence.  See  Appendix  F  for  the  complete  run  script. 
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4.0  TROUBLESHOOTING 

Table  2  provides  a  brief  troubleshooting  guide  for  when  problems  are  encountered.  It  is 
followed  by  a  more  detailed  discussion  of  each  specific  problem. 


Table  2:  GOFS  V3.0  Troubleshooting  Guide. 


Problem 

Solution 

1)  System  has  stopped 
running. 

Check  log  files  in 

/scr/ooc/ data/hycom/ GLBaO  .08/ exp  t_72 . 7/logs 

to  find  the  part  of  the  system  where  the  problem  occurred. 

2)  HYCOM  stopped  running 
due  to  negative  layer 
thicknesses,  i.e.  “neg.  dp”  is 
seen  in  the  log  file. 

a)  Reduce  the  baroclinic  time  step  (baclin)  in 
/u/home/ooc/models/hycom/GLBa0.08/expt_72.7/blkdat. input. 
If  unsuccessful,  try, 

b)  Turn  off  the  incremental  updating  (in  blkdat.input)  by 
setting  ’incflg'  =  incremental  update  flag  (0=no,  l=yes,  2=full- 
velocity). 

3)  Error  message:  Cannot  load 
ntbl  windows  on  all  selected 
nodes... 

This  is  a  machine  problem  and  resubmitting  the  script  usually 
works. 

4)  Cholesky  decomposition 
failed  error  in  a  NCODA 
subregion. 

Adjust  the  diagonal  of  the  covariance  matrices  by  changing 
the  oanl  namelist  variable  err_plus()  =  1 .  Set  it  to  a  negative 
value  after  the  analysis  has  gotten  past  the  problem,  see 
oanl.h,  the  NCODA  namelist  file. 

4.1  The  System  Has  Stopped  Running 

Explanation:  One  of  the  jobs  has  died  for  some  reason. 

Solution:  See  /scr/ooc/data/hycom/GLBa0.08/expt_72. 7/logs/. 

At  the  system  prompt,  enter 

Is  -lat  |  more 

to  see  where  the  run  died.  Check  the  appropriate  log  files  for  errors. 


4.2  GOFS  Runs  Fine  For  A  While  And  Then  Suddenly  Negative  Layer 
Thicknesses  Appear 


The  user  may  see  “neg.  dp”  in  the  log  file. 
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Explanation:  Model  is  unstable. 

Solution  1:  Reduce  the  baroclinic  time  step  (baclin)  in 

/u/home/ ooc/models/hycom/ GLBaO . 08 / expt_72 . 7 /blkdat . input 

Use  the  table  in  blkdat.input,  as  illustrated  in  Table  3  below,  to  change  the  baroclinic  time 
step  and  incstp  so  that  the  increments  are  put  in  over  6  hours. 


Table  3:  Baroclinic  time  step  and  increment  step  options  found  in  blkdat.input. 


baclin 

incstp 

75 

288 

120 

180 

150 

144 

240 

90 

300 

72 

Solution  2:  If  lowering  the  time  step  does  not  work,  the  observations  might  be  bad.  Turn 
the  assimilation  off  for  the  day.  In 

/u/home/ ooc/models/hycom/ GLBaO . 08 / expt_72 . 7 /blkdat . input 

change 

2  '  incf  lg  '  =  incremental  update  flag  (0=no,  l=yes,  2=full-velocity) 

to 

0  '  incf  lg'  =  incremental  update  flag  (0=no,  l=yes,  2=full-velocity). 

It  is  possible  that  the  user  must  go  back  to  the  previous  day,  since  the  instability  could 
have  begun  earlier. 

NOTE:  Remember  to  change  incflg  back  to  2  after  the  problem  day  has  run 
successfully, 

4.3  The  NCODA  Analysis  or  The  Model  Log  Files  Contains  An  Error 

The  error  could  be,  for  example, 

Sat  Nov  29  11:32:09  GMT  2008  [poejob]:  Cannot  load 

ntbl_windows  on  all  selected  nodes  . . . 

Explanation:  An  error  such  as  this  usually  indicates  a  machine  problem.  Resubmitting  the 
job  usually  solves  the  problem. 


9 


NRL/MR/7320— 09-9175 


GOFS  V3.0  User’s  Manual 


Solution:  Resubmit  the  script  that  failed.  If  it  is  one  of  the  NCODA  analysis  scripts, 
uncomment  the  line 

# /u /home /wallcraf /bin /msub 

$ { E } ncoda2arch_$ { idtgtod}_$ { curr_dtg } .csh  1 

in  e.g.,  727ncoda_ARCatl_2008112918_2008112518.csh  at  the  bottom  of  the  script.  This 
is  done  so  that  the  NCODA-to-HYCOM  script  is  submitted  and  the  model  run  continues. 
Make  sure  that  the  ncoda2arch  script  is  not  already  running.  If  so,  kill  the  job  before 
resubmitting  the  NCODA  analysis  for  the  subregion. 

4.4  Cholesky  Decomposition  Failed 

Explanation:  Either  the  background  error  is  too  large  or  the  observation  error  is  too 
small.  This  results  in  very  small  diagonal  elements.  With  converging  grid  locations  and 
super-obs  based  on  grid  i,j  indices,  the  data  can  be  very  close  and  the  off-diagonal 
correlations  can  be  large.  The  result  is  a  poorly  conditioned  matrix. 

Solution:  Adjust  the  diagonal  of  the  covariance  matrices  by  changing  the  oanl  namelist 
variable  err_plus()  =  1.  Set  it  to  a  negative  value  after  the  analysis  has  moved  past  the 
problem.  See  the  NCODA  namelist  file  oanl.h. 
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5.0  GOFS  V3.0  OPERATIONAL  FLOWCHART 


f  1.  Cron  job  which  starts  the  entire 
system  on  the  operational  platform 

/u/home/ooc/models/hycom/GLBa0.08/expt_90.6/mvoi/ 

906shellsub  submit  ncoda.csh 


3.  Extract  NCODA  z-levels  from 
HYCOM  archive  file 
(for  each  subregion) 

/scr/ooc/data/hycom/GLBa0.08/expt_72.7/logs/7 
27extr_${subreg}_${idtgtod}  1 8_$  {idtg}  18.com 


4.  Run  NCODA  analysis  for  each 
subregion  (run  in  parallel) 

/scr/ooc/data/hycom/GLBaO .  0  8/expt_72 . 7/logs/ 
727nocda_${subreg}_${idtgtod}  18_$  {idtg}  18. 
com 


5.  Convert  NCODA  analysis  on  z- 
levels  to  HYCOM  vertical  coordinates 

/scr/ooc/data/hycom/GLBaO .  08/expt_72 . 7/logs/ 
906ncoda2arch_${idtgtod}  18_$  {idtg}  18.com 
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6.  Run  HYCOM 

/scr/ooc/data/hycom/GLBa0.08/expt_72. 7/logs/ 
727lsf_ncoda_  ${idtgtod}  18_$  {idtg}  18.com 


7.  Make  scripts  and  start  the  next  day\ 
of  the  daily  run 

/u/home/ ooc/ models/hycom/  GLB  aO .  0  8/ expt_7  2 . 7/mvoi/ 
727shellsub_daily.csh 

Post-processing  begins  for  day  that  just  finished.  / 


E.J. Metzger  and  O.M.  Smedstad 
Naval  Research  Laboratory,  Code  7323 
January  2009 


Figure  2:  Flowchart  of  the  GOFS  Y3.0  runstream.  All  scripts  are  provided  in 
Appendices  A-F. 
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6.0  FUNCTIONAL  DESCRIPTION 

For  a  discussion  of  the  HYCOM  2.2  model  theory  and  code  description  see  the 
accompanying  Software  Design  Description  (Wallcraft  et  al.,  2009). 

7.0  NOTES 

7.1  Acronyms  and  Abbreviations 


Acronym 

Definition 

GOFS 

Global  Ocean  Forecast  System 

HYCOM 

HYbrid  Coordinate  Ocean  Model 

MPI 

Message  Passing  Interface 

NCAR 

National  Center  for  Atmospheric  Research 

NCODA 

Navy  Coupled  Ocean  Data  Assimilation 

NRL 

Naval  Research  Laboratory 

SDD 

Software  Design  Description 

UM 

User’s  Manual 

VTR 

Validation  Test  Report 
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APPENDIX  A 

Script  727shellsub_submit_ncoda.csh 

This  is  the  main  driver  script. 

#!  /bin/csh  -f 

# 

#  -  This  script  generates  the  first  set  of  scripts  needed  for  todays 

#  -  HYCOM-NCODA  run.  Most  of  the  work  is  simply  script  creation,  but 

#  -  the  last  thing  it  does  is  submit  a  job  to  create  the  forcing  needed 

#  -  by  the  ocean  model.  A  subsequent  sequence  of  scripts  gets  executed 

#  - in  turn. 

# 

#  -  Originally  scripted  by  O.M.  Smedstad,  QinetIQ-PSI,  Inc. 

#  -  ole.smedstad.ctr@nrlssc.navy.mil 

#  -  office  phone:  228-688-4365 

# 

set  path=  (  $path  /usr/bin/Xll  /u/home/wallcraf /bin  /u/home/$ { user } /bin  .  ) 
echo  $path 
# 

#  -  The  environmental  variables  ${E}  and  ${EXPT}  define  the  experiment 

#  -  number.  If  setting  up  for  a  new  experiment,  these  need  to  be  changed 

#  -  throughout  all  scripts. 

# 

setenv  E  727 

setenv  REG1  GLBa0.08 

setenv  EXPT  expt_72.7 

setenv  MVOI  /u/home/$ {user } /hycom/$ {REG1 }/$ {EXPT} /mvoi 
setenv  LOGS  /scr/$ { user } /hycom/ ${ REG1 }/${ EXPT } /logs 
# 

#  -  Make  directories  in  case  they  do  not  exist 

# 

mkdir  -p  ${L0GS} 

mkdir  -p  /scr/$ {user } /hycom/$ {REG1 }/$ {EXPT} /data/wind 
mkdir  -p  /scr/$ { user } /hycom/$ {REG1 }/${ EXPT } /data/flux 
mkdir  -p  /scr/$ { user } /hycom/$ {REG1 }/$ {EXPT } /data/ssta 
mkdir  -p  /scr/$ {user } /hycom/$ {REG1 }/$ {EXPT} /data/pcip 
mkdir  -p  /scr/$ {user } /hycom/$ {REG1 }/$ {EXPT} /data/wspd 
mkdir  -p  /scr/$ { user } /hycom/$ {REG1 }/$ {EXPT } /data/incup 
# 

#  -  The  time  of  the  NCODA  analysis 

# 

setenv  HR  18 
setenv  HR 2  ${HR} 

# 

#  -  The  day  model  is  submitted  (run)  (format  20070101) 

#  -  idtgtod:  date-time-group  for  today 

# 

setenv  idtgtod  'date  +%Y%m%d' 
fsetenv  idtgtod  $1 
# 

#  -  Go  back  5  days  from  today  to  do  a  5  day  hindcast 

#  -  idtg:  date-time-group  for  today  minus  5  days 

# 

set  idtg= ' /u/home/$ {user } /bin/addndays  yyyymmdd  ${idtgtod}  -5' 

# 
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#  -  Total  number  of  days  to  create  forcing  files 

# 

#set  maxnumdays=15 
set  maxnumdays=3 
# 

#  -  Number  of  days  of  model  integration  between  each  NCODA  analysis 

# 

setenv  nmdays  1 

# 

#  -  idtgml:  one  day  prior  to  the  system  start  date 

#  -  idtgtodml :  one  day  prior  to  the  system  start  date 

# 

setenv  idtgml  ' /u/home/$ { user } /bin/addndays  yyyymmdd  ${idtg}  -1' 
setenv  idtgtodml  ' /u/home/$ {user } /bin/addndays  yyyymmdd  ${idtgtod}  -1' 

# 

echo  'idtgtod='  ${idtgtod}  '  idtg='  ${idtg} 

# 

#  -  Start  job  at  this  time 

# 

setenv  month  'echo  ${idtgtod}  | cut  -c5-6' 
setenv  day  'echo  ${idtgtod}  | cut  -c7-8' 
setenv  nxtd  20:00 
# 

#  -  Delete  forcing  files  that  may  exist  with  idtg 

# 

/bin/rm  -f  /scr/$ {user} /hycom/$ {REG1 } /$ {EXPT} /data/ ssta/*$ { idtg} * 

/bin/rm  -f  /scr/$ {user} /hycom/$ {REG1 } /$ {EXPT} /data/ wind/ *  $ { idtg} * 

/bin/rm  -f  /scr/${user} /hycom/ $ { REG1 } / $ {EXPT} /data/ f lux/*$ { idtg} * 

/bin/rm  -f  /scr/${user} /hycom/ $ {REG1 } /$ {EXPT} /data/wspd/* $ { idtg} * 

/bin/rm  -f  /scr/$ {user} /hycom/ $ {REG1 } /$ {EXPT} /data/pcip/* $ { idtg} * 

# 

cd  /u/home/ $ { user } /hycom/ $ { REG1 } / $ { EXPT } 

# 

#  -  Make  necessary  scripts  for  each  NCODA  analysis  subregion 

#  -  for  the  given  date 

# 

foreach  REG (  ANTarc  MERatl  MERinl  MERpal  ARCatl  ARCpac  ARCocn  ) 

# 

#  -  Generate  the  script  to  plot  NCODA  fields 

# 

/bin/rm  $ {LOGS } /$ {E }plot_$ { REG}_$ { idtgtod} $ {HR}_$ { idtg} $ { HR} * . { csh, log} 
awk  -f  $ {MVOI } /MVOI . awk  tod=$ { idtgtod}  hr=${HR}  tl=${idtg}  tr=${REG}  \ 

$ { MVOI }/${E}plot. csh  >  \ 

$ { LOGS } / $ { E } plot_$ {REG }_$ { idtgtod} $ { HR}_$ { idtg } $ { HR} .csh 

# 

#  -  Generate  the  script  to  rep  NCODA  plots 

# 

/bin/rm  $ {LOGS } /$ {E }plotrcp_$ {REG}_$ { idtgtod} $ { HR}_$ { idtg} $ {HR} * . { csh, log} 
awk  -f  $ {MVOI } /MVOI . awk  hr=${HR}  tod=$ { idtgtod}  tl=${idtg}  tr=${REG}  \ 

$ {MVOI }/$ {E }plotrcp . csh  >  \ 

$ { LOGS } /$ {E }plotrcp_$ { REG }_$ { idtgtod} $ { HR}_$ { idtg } $ { HR} .csh 

# 

#  -  Generate  the  script  to  perform  the  NCODA  analysis 

# 

/bin/rm  $ {LOGS } / $ {E }ncoda_$ {REG}_$ { idtgtod} $ { HR } _ $ { idtg} $ {HR} * . { csh, log} 

awk  -f  ${ MVOI } /MVOI .awk  reg=${REG}  tod=$ { idtgtod} ${ HR}  tl=$ { idtg} $ { HR}  \ 
t2=l  t3=${idtgml }${HR}  t4=$ { E } ncoda_$ { REG}_$ { idtg} .csh  \ 

$ { MVOI } / $ { E } ncoda_$ { REG } . j  ob  >  \ 
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$ { LOGS } / $ { E } ncoda_$ { REG}_$ { idtgtod} $  { HR}_$ { idtg} $ { HR} . csh 

# 

#  Generate  the  script  to  extract  HYCOM  fields  for  NCODA  analysis 

#  todml=idtgtodml  for  the  first  analysis 

# 

/bin/rm  $ { LOGS } /$ { E } extr_$ { REG}_$ { idtgtod} $ { HR}_$ {idtg}${HR}*.{csh,log} 
awk  -f  $ {MVOI } /MVOI_daily . awk  nxtd=${nxtd}  calnmb=l  hr=${HR}  hr2=${HR} 
tl=${idtg}  \ 

tod=$ { idtgtod}  todl=$ { idtgtod}  todml=$ { idtgtodml } 

$  {MVOI } /$ { E } extr_$ {REG} . csh_zi  >  \ 

$ { LOGS } / $ { E } extr_$ {REG}_$ { idtgtod} $ { HR}_$ { idtg } $ { HR} . csh 

# 

#  Change  file  permissions 

# 

chmod  u+rx  ${ LOGS }/*${ idtg} * 

# 

end 

# 

#  -  Generate  script  to  interpolate  to  uniform  1/12  deg  grid  (GLBu0.083) 

#  and  write  netcdf  files 

# 

# /bin/rm  $ { LOGS } / $ { E } subreg_GLBuO . 083_$ { idtgtod} $ { HR } _ $ { idtg} $ { HR} * . {csh, log} 

#awk  -f  $ {MVOI } /MVOI . awk  hr=${HR}  tod=$ { idtgtod} 00  tl=$ { idtgtod} 00  \ 

#t2=l  t3=$ { idtgml } 00  $ {MVOI }/../.. /subregion/$ { E } subreg_GLBuO . 083 . com  >  \ 

#  $ { LOGS } /$ { E } subreg_GLBu0 . 083_$ { idtgtod} $ { HR}_$ { idtg } $ { HR} . csh 

# 

#  Generate  the  NCODA  prep  script  (get  all  data  for  the  period  idtg  to 

#  idtgtod)  not  necessary  for  a  real  time  job  if  the  data  are  already  on 

#  the  machine 

# 

/bin/rm  $ {LOGS } /$ { E } ncoda_PREP_$ { idtgtod} $ { HR}_$ { idtg} $ {HR} * . { csh, log} 
awk  -f  $ {MVOI } /MVOI . awk  nxtd=${nxtd}  reg=PREPl  hr=${HR}  tod=$ { idtgtod} 00  \ 
tl=$ { idtgtod} 00  t2=l  t3=$ { idtgml } 00  t4=$ { E } ncoda_$ {REG1 }_${ idtg} . csh  \ 

$ { MVOI } / $ { E } ncoda_PREP . j  ob  >  \ 

$ { LOGS } / $ { E } ncoda_PREP_$ { idtgtod } $ { HR } _$ { idtg } $ { HR } . csh 

# 

#  -  Generate  the  script  to  convert  a  HYCOM  restart  file  to 

#  a  HYCOM  archive  file 

# 

# /bin/rm  $ { LOGS } / $ { E } rest2archv_$ { idtgtod} $ { HR } _ $ { idtg} $ { HR} . csh 

#awk  -f  $ {MVOI } /MVOI . awk  nxtd=${nxtd}  hr=${HR}  tod=$ { idtgtod}  todl=$ { idtgtod}  \ 

#  tl=${idtg}  $ {MVOI } /$ { E } restart_archv_l . csh  >  \ 

#  $ {LOGS } /$ { E } rest2archv_$ { idtgtod} $ { HR } _ $ { idtg} $ {HR} .csh 

# 

#  -  Generate  the  script  to  convert  an  NCODA  file  to  a  HYCOM  archive  file 

# 

/bin/rm  $ {LOGS } /$ { E } ncoda2arch_$ { idtgtod} $ { HR}_$ { idtg} $ {HR} * . { csh, log} 

awk  -f  $ {MVOI } /MVOI . awk  hr=${HR}  tod=$ { idtgtod}  todl=$ { idtgtod}  tl=${idtg}  \ 

todml=$ { idtgtodml }  $ {MVOI }/${ E } ncoda2arch_all . csh  > 

\ 

$ { LOGS } / $ { E } ncoda2arch_$ { idtgtod} $ {HR}_$ { idtg} $ { HR} . csh 

# 

#  -  Generate  the  script  to  copy  log  files  to  newton 

# 

/bin/rm  $ {LOGS } /$ { E } cplogs_$ { idtgtod} $ { HR}_$ { idtgml } $ {HR} * . { csh, log} 
awk  -f  $ {MVOI } /MVOI . awk  hr=${HR}  tl=${ idtgml}  tod=$ { idtgtod}  \ 

$ {MVOI } / . . / cplogs . csh  >  $ { LOGS } /$ { E} cplogs_$ { idtgtod} $ {HR}_$ { idtgml } $ { HR} . csh 

# 
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#  -  Generate  the  script  to  copy  gmeta  files  to  NRL  anonymous  ftp  server 

# 

/bin/rm  $ {LOGS } /$ {E } cpgmeta_$ { idtgtod} $ {HR}_$ { idtgml } $ { HR} * . { csh, log} 
awk  -f  $ {MVOI } /MVOI . awk  hr=${HR}  tl=${ idtgml}  tod=$ { idtgtod}  \ 

$ {MVOI } / . . / cpgmeta .csh  >  $ {LOGS } /$ {E } cpgmeta_$ { idtgtod} $ { HR } $ { idtgml } $ { HR} . csh 

# 

#  -  Generate  the  script  to  copy  archv_l  files  to  newton 

# 

/bin/rm  $ {LOGS } /$ {E } cparchv_l_$ { idtgtod} $ { HR } _ $ { idtgml } $ { HR} * . { csh,  log} 

awk  -f  $ {MVOI } /MVOI . awk  hr=${HR}  tl=${ idtgml}  tod=$ { idtgtod}  \ 

$ {MVOI }/.. /cparchv_l . csh  >  \ 

$ { LOGS } /$ { E} cparchv_l_$ { idtgtod} $ { HR } _ $ { idtgml } $ {HR} .csh 

# 

#  -  Generate  the  script  to  copy  files  to  NRL  anonymous  ftp  server 

# 

/bin/rm  ${ LOGS } /ztar_$ { idtgtod} ${ HR}_$ { idtg }${ HR} .csh 

awk  -f  $ {MVOI } /MVOI . awk  hr=${HR}  hr=${HR}  tl=${idtg}  tod=$ { idtgtod}  \ 

$ {MVOI } / . . /$ { E }maktar_glbaO . 08 . com  > 

$ {LOGS } / ztar_$ { idtgtod} $ { HR } _ $ { idtg} $ {HR} .csh 

# 

#  -  Generate  the  script  to  get  the  BOGUS  files  from  file  server 

# 

/bin/rm  $ {LOGS } /getfro_$ { idtgtod} $ { HR} . sh 

awk  -f  $ {MVOI } /MVOI . awk  hr=${HR}  hr=${HR}  tl=${idtg}  tod=$ { idtgtod}  \ 

$ {MVOI } / . . /getfro . sh  >  ${ LOGS } /getfro_$ { idtgtod} ${ HR} . sh 

# 

cd  $ { LOGS } 

# 

#  -  Submit  job  getting  the  BOGUS  files 

# 

/u/home/wallcraf /bin/q_navo  getfro_$ { idtgtod} $ { HR} . sh 

# 

#  -  Generate  the  script  to  create  forcing  files  for  the  number  of 

#  -  days  in  todays  run  (maxnumdays) 

# 

/bin/rm  $ {LOGS } /$ {E } cpforcing_$ { idtgtod} $ {HR}_$ { idtg} $ { HR} * . { sh, log } 
awk  -f  $ {MVOI }/.. /HYCOM . awk  nmdays=$ {maxnumdays }  t0=$ { E }pbs_$ { idtg}  hr=${HR}  \ 
tod=$ { idtgtod}  tl=${idtg}  $ {MVOI }/../$ {E } cpforcingtomsrc_real . sh>  \ 

$ { LOGS } /$ { E} cpforcing_$ { idtgtod} $ { HR}_$ { idtg} $ {HR} . sh 

# 

/bin/rm  $ {LOGS } / $ { E }pbs_wind_prep_$ { idtgtod} $ {HR}_$ { idtg} $ { HR} * . { com,  log} 
awk  -f  $ {MVOI }/.. /HYCOM . awk  nmdays=$ {maxnumdays }  t0=$ { E }pbs_$ { idtg}  hr=${HR}  \ 
tod=$ { idtgtod}  tl=${idtg}  $ {MVOI }/../$ {E }pbs_wind_prep . com  >  \ 

$ { LOGS } / $ { E}pbs_wind_prep_$ { idtgtod} $ { HR}_$ { idtg} $ {HR} . com 

# 

#  -  Generate  the  HYCOM  run  script  for  today 

# 

/bin/rm  $ {LOGS } / $ { E }pbs_ncoda_$ { idtgtod} $ {HR}_$ { idtg} $ { HR} * . { com, log} 
awk  -f  $ {MVOI }/.. /HYCOM . awk  nmdays=$ { nmdays }  t0=$ { E }pbs_$ { idtg}  hr=${HR}  \ 
hr2=${HR}  tod=$ { idtgtod}  tl=${idtg}  $ {MVOI }/../${ E }pbs_ncoda . com  >  \ 

$ { LOGS } / $ { E}pbs_ncoda_$ { idtgtod} $ { HR}_$ { idtg} $ {HR} . com 

# 

#  -  Submit  the  forcing  preparation  job 

# 

#  on  XT 5 

/u/home/wallcraf /bin/ q_navo  $ {E } cpforcing_$ { idtgtod} $ { HR } _ $ { idtg} $ { HR} . sh 

#  on  IBM 

# /u/home/wallcraf /bin/q_navo  $ { E } pbs_wind_prep_$ { idtgtod} $ { HR}_$ { idtg} $ {HR} . com 
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# 

#  -  Submit  job  to  transfer  observations  files  for  NCODA 

# 

/u/home/wallcraf /bin/ q_navo  $ {E } ncoda_PREP_$ { idtgtod} $ { HR}_$ { idtg} $ {HR} . csh 
exit 
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APPENDIX  B 

Script  7271sf_wind_pr ep_$ {idtgtod}  1 8  $ {idtg}  18.com 

This  is  the  atmospheric  forcing  generation  script. 

#!/bin/csh  -f 
#PBS  -N  999pbs 
#PBS  -j  oe 

#PBS  -o  999pbsXX.log 
#PBS  -1  mppwidth=l 
#PBS  -1  mppnppn=l 
#PBS  -1  walltime=3 : 00 : 00 
#PBS  -W  umask=027 
#PBS  -A  NRLSS018 
#PBS  -q  internal3d 
# 

#BSUB  -J  727ncodawind 
#BSUB  -n  1 

#BSUB  -R  "span [ptile=l] " 

#BSUB  -W  4:00 
#BSUB  -q  share 
#BSUB  -P  NAVOSHYC 
# 

#  -  Job  to  create  the  forcing  needed  for  a  HYCOM-NCODA  run 

# 

set  echo 

set  timestamp 

C 

C  -  Preamble. 

C 

setenv  OS  'uname' 
switch  ($OS) 
case  1  Linux  1 : 
which  yod 

if  (!  $status)  then 
setenv  OS  XT3 
setenv  TMPDIR  / tmp 
else 

setenv  TMPDIR  /tmp 
endif 

which  aprun 
if  (!  $status)  then 
#  setenv  OS  XT4 

setenv  OS  XT5 
setenv  TMPDIR  /scr 
endif 
breaksw 
case  1 AIX ' : 
hostname 

setenv  TMPDIR  /scr/${user} 
breaksw 
default : 

echo  'Unknown  Operating  System:  '  $OS 
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echo  'configured  for  AIX  only' 
exit  (1) 
endsw 
C 

if  ( $  ?  J OBNAME )  then 

setenv  PBS_JOBNAME  $ { JOBNAME } 
setenv  PBS_JOBID  $$ 
endif 

echo  PBS_ JOBNAME  $PBS_JOBNAME  PBS_JOBID  $PBS_JOBID 
C 

C  -  E  is  expt,  P  is  permenant  directory,  S  is  /trap  directory. 

C 

# 

setenv  E  727 
setenv  EXPT  expt_72.7 
setenv  REG  GLBa0.08 
# 

#  -  These  variables  are  set  by  the  awk  script 

# 

setenv  nmdays  10 
setenv  idtg  20070209 

setenv  idtg2  ' /u/home/$ {user } /bin/addndays  yyyymmdd  ${idtg}  +1' 

setenv  idtgtod  20070214 

setenv  HR  00 

setenv  HR2  00 

# 

echo  ${idtg}  ${idtgtod} 

# 

set  LOGS=/scr/$ {user} /hycom/ $ {REG} /$ {EXPT } /logs 
mkdir  -p  ${LOGS} 

set  RUN=/ u/home/ $ {user } /hycom/ $ {REG} /$ { EXPT } 
cd  $ { RUN } 

# 

setenv  P  $cwd 
switch  ($OS) 
case  'AIX': 
case  ' XT5 ' : 

# 

#  -  substitute  /scr  for  /u/home 

# 

setenv  S  'echo  $cwd  |  awk  '{print  "/scr"  substr ($0, 8, length) }' ' 
breaksw 
case  1 XT3 ’ : 
case  ' XT4 ' : 
case  ' OSF1 '  : 

#  substitute  /work  for  /??????? 

setenv  S  'echo  $cwd  |  awk  '{print  "/work"  substr  ( $0 , 3 , length)  }' 
breaksw 
case  ' IRIX64 ' : 

#  substitute  /scr  for  /u/home 

setenv  S  'echo  $cwd  |  awk  '{print  "/scr"  substr ($0, 8, length) }' ' 
breaksw 
case  ' unicos ' : 

#  substitute  /trap  for  /u/b 

setenv  S  'echo  $cwd  |  awk  '{print  "/trap"  substr ($0, 5, length) }' ' 
breaksw 


endsw 

C 
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Is  -laFq 

set  idtgtodl= ' echo  ${idtgtod}  |  cut  -cl-8' 

C 

C  -  check  the  RUNNING  flag. 

C 

if  (  -e  RUNNING  &&  !  -e  RUNNING_$PBS_JOBID)  then 
C 

C  -  MODEL  IS  ALREADY  RUNNING  -  EXIT. 

C 

exit 

endif 

touch  RUNNING 

touch  RUNNING_$PBS_JOBID 

C 

C  -  Generate  the  next  model  script. 

C 

setenv  Y01  103 
setenv  AB  a 

# 

setenv  SCRIPT  $ { E } y$ { idtg } . com 
/bin/rm  -f  ${SCRIPT} 

awk  -f  $ { RUN} /$ {E } . awk  nmdays=$ { nmdays }  y01=${Y01}  hr=${HR}  tod=$ { idtgtod}  \ 
ab=${AB}  td=${idtg}  $ {RUN} /$ {E }_wind. com  >  ${SCRIPT} 

# 

#  -  Run  the  Script. 

# 

set  script  =  $SCRIPT 

set  reqname  =  $ { PBS_JOBNAME } 

In  -fs  ${ reqname }. log  $script : r . log 

# 

mkdir  -p  $S 

cp  $ { SCRIPT }  $S/ $ { SCRIPT } 
cd  $S 
# 

# - 

csh  $ { SCRIPT } 

# - 

# 

Is  -laFq 

# 

cd  $P 

# 

#  -  Final  Clean  Up. 

# 

/bin/rm  -f  RUNNING 

/bin/rm  -f  RUNNING_$PBS_JOBID 

# 

#  -  Make  links  for  the  forcing  files  created  today 

# 

while (${ idtg2 }  <=  ${idtgtodl}) 

# 

cd  /scr/$ {user } /hycom/ $ {REG} / $ { EXPT } /data /wind 
/bin/rm  *${idtg2}.* 

In  -s  tauewd_$ { idtg} . a  tauewd_$ { idtg2 } . a 
In  -s  tauewd_$ { idtg} . b  tauewd_$ { idtg2 } . b 
In  -s  taunwd_$ { idtg} . a  taunwd_$ { idtg2 } . a 
In  -s  taunwd_$ { idtg} . b  taunwd_$ { idtg2 }  .  b 
# 
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cd  /scr/$ {user} /hycom/$ {REG} / $ { EXPT } / data /flux 
/bin/rm  *${idtg2}.* 

In  -s  airtmp_$ { idtg} . a  airtmp_$ { idtg2 } . a 
In  -s  air tmp_$ { idtg }. b  airtmp_$ { idtg2 } . b 
In  -s  radf lx_$ { idtg } . a  radf lx_$ { idtg2 } . a 
In  -s  radf lx_$ { idtg }. b  radf lx_$ { idtg2 } . b 
In  -s  shwf lx_$ { idtg } . a  shwf lx_$ { idtg2 } . a 
In  -s  shwf lx_$ { idtg }. b  shwf lx_$ { idtg2 } . b 
In  -s  vapmix_$ { idtg } . a  vapmix_$ { idtg2 } . a 
In  -s  vapmix_$ { idtg } . b  vapmix_$ { idtg2 } . b 

cd  /scr/$ {user } /hycom/ $ {REG} /$ { EXPT } /data/ssta 
/bin/rm  *${idtg2}.* 

In  -s  surtmp_$ { idtg} . a  surtmp_$ { idtg2 } . a 
In  -s  surtmp_$ { idtg} . b  surtmp_$ { idtg2 } . b 

cd  /scr/$ {user } /hycom/ $ {REG} / $ { EXPT } / data/pcip 
/bin/rm  *${idtg2}.* 

In  -s  precip_$ { idtg} . a  precip_$ { idtg2 } . a 
In  -s  precip_$ { idtg} .b  precip_$ { idtg2 } . b 

cd  /scr/$ {user } /hycom/ $ {REG} /$ { EXPT } / data/wspd 
/bin/rm  *${idtg2}.* 

In  -s  wndspd_$ { idtg} . a  wndspd_$ { idtg2 } . a 
In  -s  wndspd_$ { idtg} . b  wndspd_$ { idtg2 } . b 
# 

set  idtg2=' /u/home/$ { user } /bin/addndays  yyyymmdd  ${idtg2}  +1' 
end 

cd  $ { LOGS } 

# 

#  Submit  the  job  extracting  an  archive  file  from  a  restart  file. 

#  The  job  will  wait  for  the  time  slot  dedicated  to  the  HYCOM  run 

# 

#/u/home/wallcraf /bin/q_navo  $ { E} rest2archv_$ { idtgtod} $ {HR}_$ { idtg} $ { HR} . csh 

/u/home/wallcraf /bin/ q_navo  $ {E } ncoda_PREP_$ { idtgtod} $ { HR } _ $ { idtg} $ {HR} .csh 

# 

#  - Exit. 

# 

exit 
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APPENDIX  C 

Script  727extr_${subreg}_${idtgtod}18_${idtg}18.com 

This  script  extracts  NCODA  z-levels  from  the  HYCOM  archive  file  for  each  NCODA 
subregion. 


#!/bin/csh  -f 
#PBS  -N  999pbs 
#PBS  -j  oe 

#PBS  -o  999pbsXX.log 
#PBS  -1  mppwidth=l 
#PBS  -1  mppnppn=l 
#PBS  -1  walltime=l : 00 : 00 
#PBS  -W  umask=027 
#PBS  -A  NRLSS018 
#PBS  -q  internal3d 
# 

#BSUB  -J  pmast 

#BSUB  -n  1 

#BSUB  -W  0:30 

#BSUB  -R  "span [ptile=l] " 

#BSUB  -q  challenge 
#BSUB  -P  NRLSSC3 J 
#####BSUB  -q  internal3d 
####BSUB  -b  10:24:20:00 
# 

date 

set  echo 

set  pget=~wallcraf /bin/pget 

# 

set  OS='uname' 
switch  ($OS) 
case  1  Linux  1 : 

which  aprun 

if  (!  $status)  then 

set  APRUN=' aprun  -n  1  ' 
set  SRC=~wallcraf /hycom/ALLcnl 
else 

set  APRUN= ' ' 

set  SRC=~wallcraf /hycom/ALL 
endif 
breaksw 
case  ' AIX ' 
set  APRUN= ' ' 

set  SRC=~wallcraf /hycom/ALL 
breaksw 
default : 
set  APRUN= ' ' 

set  SRC=~wallcraf /hycom/ALL 
endsw 
# 

#  -  extract  3-d  fields  from  a  single  HYCOM  archive  file. 

# 

setenv  CALNMB  1 
setenv  HR  18 
setenv  HR2  18 
setenv  idtg  20080501 
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setenv  idtgtod  20080501 
setenv  idtgtodpl  20080501 
setenv  idtgtodml  20080501 

set  idtgout=' /u/home/$ { user } /bin/addndays  yyyymmdd  ${idtg}  -1' 

# 

echo  'DATES  TO  EXTRACT  '  ${idtg} 

# 

set  EXPTN=expt_72 . 7 
set  REGN=MERpal 0 . 08 
# 

setenv  E  727 
set  EXPT=expt_72 . 7 
set  REG=GLBa0708 
set  CASE=nowcast 
setenv  T  09 
# 

set  maxdepth=002500 
set  surface=000000 
set  forecast=00240000 
# 

set  TOPO=/u/home/$ { user } /hycom/ $ { REG} /topo 

set  SCR=/scr/$ { user } /hycom/ ${ REG} /$ {EXPT} 

set  LOGS=/scr/$ {user} /hycom/ $ {REG} /$ { EXPT } /logs 

set  SCRN=/scr/ $ {user} /hycom/ $ {REGN} / $ {EXPTN } 

mkdir  -p  /scr/$ {user } /hycom/$ {REG} /topo 

mkdir  -p  ${SCR}  ${SCRN} 

set  INP=$ { SCR} / $ { CASE } 

set  OUT=$ { SCRN} / raw/ $ { CASE } 

mkdir  -p  ${OUT} 

echo  $ { SCR}  $ { SCRN} 

# 

cd  $ { OUT } 

/bin/rm  ${OUT}/*[ab] 

# 

goto  TOPO 

touch  $ {TOPO} /regional . depth. a  ${ TOPO } /regional . depth . b 
if  ( - z  ${TOPO}/regional .depth. a)  then 

${pget}  ${ D} /../.. /topo/depth_$ { REG}_$ {T }. a  ${ TOPO} /regional . depth . a  & 
endif 

if  ( - z  ${ TOPO } /regional . depth . b)  then 

${pget}  ${ D} /../.. /topo/depth_$ { REG}_$ {T }. b  ${TOPO}/regional .depth. b  & 
endif 

touch  ${ TOPO} /regional . grid . a  $ {TOPO} /regional . grid. b 
if  ( - z  ${ TOPO } /regional . grid. a)  then 

${pget}  ${ D} /../.. /topo/regional . grid. a  ${TOPO}/regional .grid. a  & 
endif 

if  ( - z  ${ TOPO } /regional . grid. b)  then 

${pget}  ${ D} /../.. /topo/regional . grid. b  ${ TOPO } /regional . grid. b  & 
endif 

touch  $ {TOPO} /regional_mask.a 

if  ( —  z  $ { TOPO } /regional_mask . a)  then 

${pget}  ${ D} /../.. /topo/landsea_$ { REG} . a  $ {TOPO} /regional^mask . a  & 
endif 
wait 
TOPO: 

# 

#ln  -s  ${ TOPO} /depth_${ REG }_${T} .a  regional . depth . a 
#ln  -s  ${ TOPO} /depth_${ REG }_${T} .b  regional . depth . b 
#ln  -s  ${TOPO}/regional .grid. a  . 

#ln  -s  ${TOPO}/regional .grid.b  . 

#ln  -s  ${TOPO}/regional_mask.a  . 

# 

In  -s  ${ INP }/.. /data/regional . depth . a  . 
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In  -s  ${ INP} /.. /data/regional . depth .b  . 

In  -s  ${ INP} /.. /data/regional . grid. a  . 

In  -s  ${ INP} /.. /data/regional . grid. b  . 

#/bin/cp  ${TOPO}/regional_mask.a  . 

# 

#/bin/cp  $ { TOPO} /depth_$ {REG}_$ {T }  .  a  regional . depth . a 
#/bin/cp  $ { TOPO} /depth_$ {REG}_$ {T } . b  regional . depth . b 
#/bin/cp  ${ TOPO} /regional . grid. a  . 

#/bin/cp  ${ TOPO} /regional . grid. b  . 

#/bin/cp  ${TOPO}/regional_mask.a  . 

# 

setenv  FOR022A  ${ OUT } /hycom2dath_000000_$ { idtgout }${ HR}_$ { forecast } .a 
setenv  FOR022  ${ OUT } /hycom2dath_000000_$ { idtgout }${ HR}_$ { forecast } .b 
setenv  FOR033A  ${ OUT } /hycom3dt_$ {maxdepth}_$ { idtgout }${ HR}_$ { forecast } .a 
setenv  FOR033  $ { OUT} /hycom3dt_$ {maxdepth}_$ { idtgout} ${ HR }_${ forecast } . b 
setenv  FOR034A  ${ OUT } /hycom3ds_$ {maxdepth}_$ { idtgout }${ HR}_$ { forecast } .a 
setenv  FOR034  $ { OUT} /hycom3ds_$ {maxdepth }_$ { idtgout} $ { HR }_$ { forecast } . b 

#  setenv  FOR035A  ${ OUT } /hycom3dden_$ {maxdepth }_${ idtgout }${ HR}_$ { forecast }. a 

#  setenv  FOR035  ${ OUT } /hycom3dden_$ {maxdepth}_$ { idtgout }${ HR}_$ { forecast }. b 
setenv  FOR037A  ${ OUT } /hycom3du_$ {maxdepth}_$ { idtgout }${ HR}_$ { forecast } .a 
setenv  FOR037  $ { OUT} /hycom3du_$ {maxdepth }_$ { idtgout} $ { HR }_$ { forecast } . b 
setenv  FOR038A  ${ OUT } /hycom3dv_$ {maxdepth}_$ { idtgout }${ HR}_$ { forecast } .a 
setenv  FOR038  $ { OUT} /hycom3dv_$ {maxdepth}_$ { idtgout} $ { HR }_$ { forecast } . b 

#  setenv  FOR040A  ${ OUT } /hycom3dp_$ {maxdepth }_${ idtgout }${ HR}_$ { forecast }. a 

#  setenv  FOR040  ${ OUT } /hycom3dp_$ {maxdepth}_$ { idtgout }${ HR}_$ { forecast }. b 
/bin/rm  $FOR033A  $FOR034A  $FOR037A  $FOR038A 

/bin/rm  $FOR033  $FOR034  $FOR037  $FOR038 
/bin/rm  $FOR022 
/bin/rm  $FOR022A 

# 

${APRUN}  $ { SRC } /archive/src/archv2data3z  «E-o-D 
$ { INP} / archv. $ { idtgtodml } $ { HR } _ $ { idtg} $ { HR2 } . a 


HYCOM 

000 

' iexpt  ' 

=  experiment  number  xlO  (000=from  archive  file) 

3 

' yr f lag ' 

=  days  in  year  flag  ( 0=3 60 J1 6 , 1=36 6 J1 6 , 2=3 66 J01 , 3 

4500 

'  idm 

'  =  longitudinal  array  size 

3298 

'  j  dm 

'  =  latitudinal  array  size 

32 

'  kdm 

'  =  number  of  layers 

34.0  'thbase'  =  reference  density  (sigma  units) 

0  'smooth'  =  smooth  the  layered  fields  (0=F,1=T) 

1  'baclin'  =  extract  baroclinic  velocity  ( 0=total , l=baroclinic) 

1  'xyward'  =  output  original  unrotated  velocities  ( 0=no : DEFAULT, l=yes) 

199  ' iorign '  =  i-origin  of  plotted  subregion 

360  'jorign'  =  j-origin  of  plotted  subregion 

2525  ' idmp  '  =  i-extent  of  plotted  subregion  (<=idm;  0  implies  idm) 

1841  'jdmp  '  =  j-extent  of  plotted  subregion  (<=jdm;  0  implies  jdm) 

2  ' itype  '  =  interpolation  type  ( 0=sample, l=linear, 2=parabolic) 


42 

'  kzi 

i  — 

number  of 

depths  to  sample 

0.0 

'  zi 

i 

= 

sample 

cell 

interfaces 

1 

3.5 

'  zi 

i 

= 

sample 

cell 

interfaces 

2 

6.5 

'  zi 

i 

= 

sample 

cell 

interfaces 

3 

9.5 

'  zi 

i 

= 

sample 

cell 

interfaces 

4 

14 . 5 

'  zi 

i 

= 

sample 

cell 

interfaces 

5 

19.5 

'  zi 

i 

= 

sample 

cell 

interfaces 

6 

30.5 

'  zi 

i 

= 

sample 

cell 

interfaces 

7 

41.5 

'  zi 

i 

= 

sample 

cell 

interfaces 

8 

58.5 

'  zi 

i 

= 

sample 

cell 

interfaces 

9 

87.5 

'  zi 

i 

= 

sample 

cell 

interfaces 

10 

112 . 5 

'  zi 

i 

= 

sample 

cell 

interfaces 

11 

137.5 

'  zi 

i 

= 

sample 

cell 

interfaces 

12 

162.5 

'  zi 

i 

= 

sample 

cell 

interfaces 

13 

187.5 

'  zi 

i 

= 

sample 

cell 

interfaces 

14 

212 . 5 

'  zi 

i 

= 

sample 

cell 

interfaces 

15 
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237.5 

'  zi 

'  = 

sample 

cell 

interfaces 

16 

262.5 

'  zi 

'  = 

sample 

cell 

interfaces 

17 

287.5 

'  zi 

'  = 

sample 

cell 

interfaces 

18 

312.5 

'  zi 

'  = 

sample 

cell 

interfaces 

19 

375.0 

'  zi 

'  = 

sample 

cell 

interfaces 

20 

425.0 

'  zi 

'  = 

sample 

cell 

interfaces 

21 

475.0 

'  zi 

'  = 

sample 

cell 

interfaces 

22 

525.0 

'  zi 

'  = 

sample 

cell 

interfaces 

23 

650.0 

'  zi 

'  = 

sample 

cell 

interfaces 

24 

750.0 

'  zi 

'  = 

sample 

cell 

interfaces 

25 

850.0 

'  zi 

'  = 

sample 

cell 

interfaces 

26 

950.0 

'  zi 

'  = 

sample 

cell 

interfaces 

27 

1050.0 

'  zi 

'  = 

sample 

cell 

interfaces 

28 

1150.0 

'  zi 

'  = 

sample 

cell 

interfaces 

29 

1250.0 

'  zi 

'  = 

sample 

cell 

interfaces 

30 

1350.0 

'  zi 

'  = 

sample 

cell 

interfaces 

31 

1450.0 

'  zi 

'  = 

sample 

cell 

interfaces 

32 

1550.0 

'  zi 

'  = 

sample 

cell 

interfaces 

33 

1650.0 

'  zi 

'  = 

sample 

cell 

interfaces 

34 

1750.0 

'  zi 

'  = 

sample 

cell 

interfaces 

35 

1850.0 

'  zi 

'  = 

sample 

cell 

interfaces 

36 

1950.0 

'  zi 

'  = 

sample 

cell 

interfaces 

37 

2050.0 

'  zi 

'  = 

sample 

cell 

interfaces 

38 

2150.0 

'  zi 

'  = 

sample 

cell 

interfaces 

39 

2250.0 

'  zi 

'  = 

sample 

cell 

interfaces 

40 

2350.0 

'  zi 

'  = 

sample 

cell 

interfaces 

41 

2450.0 

'  zi 

'  = 

sample 

cell 

interfaces 

42 

2550.0 

'  zi 

'  = 

sample 

cell 

interfaces 

43 

0 

'botio  ' 

=  bathymetry 

I/O 

unit  (0  no 

1/ 

22 

' athio ' 

=  average  density  I/O  unit  (0  no  I/O) 

0 

'mltio  ' 

=  mix . 1 . thk . 

I/O  unit  (0  no  I/O) 

o 

o 

' tempml 

'  =  temperature  jump  across  mixed-layer 

(degC, 

0  no  I/O) 

o 

o 

' densml 

'  =  density  jump  across  mixed-layer 

(kg/m3. 

0  no  I/O) 

0 

'infio  ' 

=  intf.  depth 

I/O  unit  (0  no  I/O,  <0  label  with 

layer  #) 

0 

' wvlio 

'  =  w-velocity 

I/O  unit  (0  no  I/O) 

37 

'uvlio  ' 

=  u-velocity 

I/O  unit  (0  no  I/O) 

38 

'vvlio  ' 

=  v-velocity 

I/O  unit  (0  no  I/O) 

0 

' splio 

'  =  speed 

I/O  unit  (0  no  I/O) 

33 

' temio 

'  =  temperature  I/O  unit  (0  no  I/O) 

34 

' salio 

'  =  salinity 

I/O  unit  (0  no  I/O) 

0 

' tthio 

'  =  density 

I/O  unit  (0  no  I/O) 

E-o-D 

# 

#/usr/lpp/LoadL/full/bin/llq  -w  $LOADL_STEP_ID 

# 

#  convert  HYCOM  .a  files  to  RAW  files  (no  padding,  spval=l.elO) . 

#  comment  this  out  if  you  don't  need  RAW  files. 

# 

foreach  t  (  3dt  3ds  3du  3dv  ) 

if  (-e  ${ OUT } /hycom$ { t}_$ {maxdepth}_$ { idtgout }$ {HR}_$ { forecast }. a)  then 
/bin/rm  -f  ${ OUT } /hycom$ { t }_$ {maxdepth}_$ { idtgout }${ HR}_$ { forecast } .A 
${APRUN}  $ { SRC } /bin/hycom2raw 

$ {OUT } /hycom$ { t }_$ {maxdepth}_$ { idtgout }${ HR}_$ { forecast }. a  2525  1841  -999.0 
$ {OUT } /hycom$ { t }_$ {maxdepth }_$ { idtgout} $ { HR }_${ forecast } .A 
endif 
end 


# 

cd  $ { OUT } 

# 

#  -  extract  2-d  fields  from  a  single  HYCOM  archive  file. 

# 

set  idtgin=$ { idtg} 

#  foreach  HR2  (  21  00  03  06  09  12  15  18  ) 
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foreach  HR2  (  18  ) 
if (  $ { HR2 }  ==  21  )  then 
set  forecast=00030000 

set  idtgin=' /u/home/ $ { user } /bin/addndays  yyyymmdd  ${idtg}  -1' 
if (  $ { CALNMB }  ==  1  )  then 

set  idtgtodin=' /u/home/$ { user } /bin/addndays  yyyymmdd  ${idtgtod}  -1' 
else 

set  idtgtodin=$ { idtgtod} 
endif 
endif 

if (  $ { HR2 }  ==  00  )  then 
set  forecast=00060000 
set  idtgin=$ { idtg} 
if (  ${ CALNMB}  ==  1  )  then 

set  idtgtodin=' /u/home/$ { user } /bin/addndays  yyyymmdd  ${idtgtod}  -1' 
else 

set  idtgtodin=$ { idtgtod} 
endif 
endif 

if (  $ { HR2 }  ==  03  )  then 
set  forecast=00090000 
set  idtgin=$ { idtg} 
if (  ${ CALNMB}  ==  1  )  then 

set  idtgtodin=' /u/home/$ { user } /bin/addndays  yyyymmdd  ${idtgtod}  -1' 
else 

set  idtgtodin=$ { idtgtod} 
endif 
endif 

if (  $ { HR2 }  ==  06  )  then 
set  forecast=00120000 
set  idtgin=$ { idtg} 
if (  ${ CALNMB}  ==  1  )  then 

set  idtgtodin=' /u/home/$ { user } /bin/addndays  yyyymmdd  ${idtgtod}  -1' 
else 

set  idtgtodin=$ { idtgtod} 
endif 
endif 

if (  $ { HR2 }  ==  09  )  then 
set  forecast=00150000 
set  idtgin=$ { idtg} 
if (  ${ CALNMB}  ==  1  )  then 

set  idtgtodin=' /u/home/$ { user } /bin/addndays  yyyymmdd  ${idtgtod}  -1' 
else 

set  idtgtodin=$ { idtgtod} 
endif 
endif 

if (  $ { HR2 }  ==  12  )  then 
set  forecast=00180000 
set  idtgin=$ { idtg} 
if (  ${ CALNMB}  ==  1  )  then 

set  idtgtodin=' /u/home/$ { user } /bin/addndays  yyyymmdd  ${idtgtod}  -1' 
else 

set  idtgtodin=$ { idtgtod} 
endif 
endif 

if (  $ { HR2 }  ==  15  )  then 
set  forecast=00210000 
set  idtgin=$ { idtg} 
if (  ${ CALNMB}  ==  1  )  then 

set  idtgtodin=' /u/home/ ${ user } /bin/addndays  yyyymmdd  ${idtgtod}  -1' 
else 

set  idtgtodin=$ { idtgtod} 
endif 
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endif 

if (  $ { HR2 }  ==  18  )  then 
set  forecast=00240000 
set  idtgin=$ { idtg} 
set  idtgtodin=$ { idtgtodml } 
endif 
# 

#  setenv  FOR022A  ${ OUT } /hycom2dSf sd_$ { surface }_${ idtgout }${ HR}_$ { forecast }. a 

#  setenv  FOR022  ${ OUT } /hycom2dSf sd_$ { surf ace }_${ idtgout }${ HR}_$ { forecast }. b 
setenv  FOR023A  $ {OUT } /hycom2df sd_$ { surface }_$ { idtgout} $ { HR} _${ forecast } .a 
setenv  FOR023  ${ OUT } /hycom2df sd_$ { surf ace }_${ idtgout }${ HR}_$ { forecast } .b 

#  setenv  FOR024A  ${ OUT } /hycom2dNf sd_$ { surface }_${ idtgout }${ HR}_$ { forecast }. a 

#  setenv  FOR024  ${ OUT } /hycom2dNf sd_$ { surf ace }_${ idtgout }${ HR}_$ { forecast }. b 
setenv  FOR025A  $ {OUT } /hycom2dt_$ { surface }_$ { idtgout} $ { HR }_${ forecast } .a 
setenv  FOR025  ${ OUT } /hycom2dt_$ { surf ace }_${ idtgout }${ HR}_$ { forecast } .b 
setenv  FOR026A  $ {OUT } /hycom2dice_$ { surface }_$ { idtgout} $ { HR }_${ forecast } .a 
setenv  FOR026  ${ OUT } /hycom2dice_$ { surf ace }_${ idtgout }${ HR}_$ { forecast } .b 
setenv  FOR027A  $ {OUT } /hycom2dmix_$ { surface }_$ { idtgout} $ { HR }_${ forecast } .a 
setenv  FOR027  ${ OUT } /hycom2dmix_$ { surf ace }_${ idtgout }${ HR}_$ { forecast } .b 
/bin/rm  $FOR023A  $FOR025A  $FOR026A  $FOR027A 

/bin/rm  $FOR023  $FOR025  $FOR026  $FOR027 

# 

${APRUN}  $ { SRC } /archive/src/archv2data2d  «E-o-D 

$ { INP} / archv. $ { idtgtodin} $ { HR } _ $ { idtgin } $ { HR2 } . a 

HYCOM 

000  ' iexpt  '  =  experiment  number  xlO  (000=from  archive  file) 


3 

' yr f lag 

1 

=  days  in  year  flag  (0=360J16 

,1= 

366J16,2=366J01,3= 

=actual ) 

4500 

'  idm 

1 

=  longitudinal  array 

size 

3298 

'  j  dm 

1 

=  latitudinal  array 

size 

1 

'  kdm 

1 

=  number  of  layers 

34 . 0 

'  thbase ' 

= 

reference  density  (sigma 

units) 

0 

' smooth ' 

= 

smooth  fields  before 

plotting  (0=F,1=T) 

0 

'mthin  ' 

= 

mask  thin  layers 

from  plots 

(0=F, 1=T ) 

1 

' xyward 

=  output  original 

unrotated  velocities  ( 0=no : DEFAULT, l=yes 

199 

' iorign 

=  i-origin  of  plotted 

subregion 

360 

'  j  orign 

=  j-origin  of  plotted 

subregion 

2525 

'  idmp 

=  i-extent  of  plotted 

subregion 

(<=idm;  0  implies 

idm) 

1841 

'  jdmp 

=  j-extent  of  plotted 

subregion 

(<=jdm;  0  implies 

j  dm) 

0 

'botio  ' 

= 

bathymetry 

I/O 

unit 

(0 

no 

I/O) 

0 

' flxio  ' 

= 

surf,  heat  flux 

I/O 

unit 

(0 

no 

I/O) 

0 

'empio  ' 

= 

surf,  evap-pcip 

I/O 

unit 

(0 

no 

I/O) 

0 

' ttrio  ' 

= 

surf,  temp  trend 

I/O 

unit 

(0 

no 

I/O) 

0 

' strio  ' 

= 

surf,  sain  trend 

I/O 

unit 

(0 

no 

I/O) 
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' icvio  ' 

= 

ice  coverage 

I/O 

unit 

(0 

no 

I/O) 

0 

' ithio  ' 

= 

ice  thickness 

I/O 

unit 

(0 

no 

I/O) 

0 

' ictio  ' 

= 

ice  temperature 

I/O 

unit 

(0 

no 

I/O) 

23 

' sshio  ' 

= 

sea  surf,  height 

I/O 

unit 

(0 

no 

I/O) 

0 

'bsfio  ' 

= 

baro.  strmfn. 

I/O 

unit 

(0 

no 

I/O) 

0 

'uvmio  ' 

= 

mix.  lay.  u-vel. 

I/O 

unit 

(0 

no 

I/O) 

0 

'vvmio  ' 

= 

mix.  lay.  v-vel. 

I/O 

unit 

(0 

no 

I/O) 

0 

' spmio  ' 

= 

mix.  lay.  speed 

I/O 

unit 

(0 

no 

I/O) 

0 

'bltio  ' 

= 

bnd.  lay.  thick. 

I/O 

unit 

(0 

no 

I/O) 
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'mltio  ' 

= 

mix.  lay.  thick. 

I/O 

unit 

(0 

no 

I/O) 

0 

'sstio  ' 

= 

mix.  lay.  temp. 

I/O 

unit 

(0 

no 

I/O) 

0 

' sssio  ' 

= 

mix.  lay.  sain. 

I/O 

unit 

(0 

no 

I/O) 

0 

'ssdio  ' 

= 

mix.  lay.  dens. 

I/O 

unit 

(0 

no 

I/O) 

-1  ' kf  '  =  layer  to  plot  (=0  end  layer  plots;  <0  label  with  layer  #) 

1  ' kl  '  =  last  output  layer 

0  'uvlio  '  =  layer  k  u-vel .  I/O  unit  (0  no  I/O) 

0  'vvlio  '  =  layer  k  v-vel .  I/O  unit  (0  no  I/O) 

0  'splio  '  =  layer  k  speed.  I/O  unit  (0  no  I/O) 

0  'infio  '  =  layer  k  i.dep.  I/O  unit  (0  no  I/O) 

0  ' thkio  '  =  layer  k  thick.  I/O  unit  (0  no  I/O) 
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25  ' temio  '  =  layer  k  temp  I/O  unit  (0  no  I/O) 

0  ' salio  '  =  layer  k  sain.  I/O  unit  (0  no  I/O) 

0  'tthio  '  =  layer  k  dens,  I/O  unit  (0  no  I/O) 

0  ' sfnio  '  =  layer  k  strmfn.  I/O  unit  (0  no  I/O) 

0  'kf  '  =  first  output  layer  (=0  end  output;  <0  label  with  layer  #) 

E-o-D 

# 

#  convert  ice  concentration  to  percent 

# 

/bin/rm  hycom2dice_$ { surf ace }_${ idtgout }${ HR}_$ { forecast } .a2 
${APRUN}  $ { SRC } /bin/hycom_expr 

hycom2dice_$ { surface}_$ { idtgout }$ {HR}_$ { forecast }. a  ONE  2525  1841  100.0  0.0 

hycom2dice_$ { surface }_$ { idtgout } $ { HR } _ $ { forecast} . a2 

/bin/mv  hycom2dice_$ { surf ace }_${ idtgout }${ HR}__$ { forecast } .a2 

hycom2dice_$ { surface }_$ { idtgout } $ { HR } _ $ { forecast} . a 

# 

foreach  t  (  fsd  t  ice  mix) 

if  (-e  ${ OUT} /hycom2d$ { t}_$ { surface }_${ idtgout }${ HR}_$ { forecast }. a)  then 
/bin/rm  -f  $ { OUT} /hycom2d$ { t }_${ surface }_$ { idtgout} $ { HR }_${ forecast } .A 
${APRUN}  $ { SRC } /bin/hycom2raw 

$ {OUT } /hycom2d$ { t }_${ surface}_$ { idtgout }$ {HR}_$ { forecast }. a  2525  1841  -999.0 

$ {OUT } /hycom2d$ { t } _ $ { surface }_$ { idtgout } $ { HR } _ $ { forecast} .A 

endif 

end 

# 

#/usr/lpp/LoadL/full/bin/llq  -w  $LOADL_STEP_ID 

# 

#  calculate  steric  ssh 

# 

set  MEAN=  0 9 4_ar  chMNA . 0011_001 5_ATH_$ { REGN } .a 

/bin/cp  / u/home/ ${user}/hycom/$ {REG } /means td/ 094_archMNA. 0011_001 5_ATH_$ { REGN } .a 

# 

foreach  a  (  000000_$ { idtgout }${ HR}_$ { forecast }  ) 

/bin/rm  -f  hycom2d [NS ] f sd* . ? 

${APRUN}  $ { SRC } /bin/hycom_ster icssh  hycom2dath_$ { a } . a  ${MEAN}  2525  1841 
hycom2dSfsd_${a} .a  >!  hycom2dSf sd_$ { a } . b 

${APRUN}  $ { SRC } /bin/hycom_expr  hycom2df sd_$ { a } . a  hycom2dSf sd_$ { a } . a  2525 

1841  1.0  -1.0  hycom2dNf sd_$ { a } . a  >!  hycom2dNf sd_$ { a } . b 
end 
# 

foreach  t  (  Nfsd  Sfsd) 

if  (-e  ${ OUT} /hycom2d$ { t}_$ { surface }_${ idtgout }${ HR}_$ { forecast }. a)  then 
/bin/rm  -f  $ { OUT} /hycom2d$ { t }_${ surface }_$ { idtgout} $ { HR }_${ forecast } .A 
${APRUN}  $ { SRC } /bin/hycom2raw 

$ {OUT } /hycom2d$ { t }_${ surface}_$ { idtgout }$ {HR}_$ { forecast }. a  2525  1841  -999.0 

$ {OUT } /hycom2d$ { t } _ $ { surface }_$ { idtgout } $ { HR } _ $ { forecast} .A 

endif 

end 

# 

#  forecash  HR2 
end 

# 

date 

# 

#  submit  ncoda  analysis 

# 

cd  $ { LOGS } 

~wallcraf /bin/q_navo  $ { E } ncoda_MERpal_$ { idtgtodpl } $ { HR } _ $ { idtg} $ { HR} . csh 
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APPENDIX  D 

Script  727nocda_${subreg}_${idtgtod}18_${idtg}18.com 

This  is  a  script  to  run  the  NCODA  analysis  for  each  subregion. 


# ! /bin/ksh 
#PBS  -N  999pbs 
#PBS  -j  oe 

#PBS  -o  999pbsXX.log 
#PBS  -1  mppwidth=416 
#PBS  -1  mppnppn=8 
#PBS  -1  walltime=4 : 00 : 00 
#PBS  -W  umask=027 
#PBS  -A  NRLSS018 
#PBS  -q  internal3d 
# 

#BSUB  -J  727ncoda_MERpal . job 
#BSUB  -n  192 
#BSUB  -a  poe 
#BSUB  -W  3:30 

#BSUB  -R  "rusage [ntbl_windows=32 ] span [ptile=l 6 ] " 

#BSUB  -q  challenge 
#BSUB  -P  NRLSSC3 J 
#####BSUB  -q  internal3d 
#####BSUB  -P  NAVOSHYC 
#BSUB  -e  err_MERpal . % J 
# 

#  -  set  days  to  run 

# 

inumd=l 

# 

#  -  analysis  date 

# 

idtg=2008050118 

# 

#  -  forecast  from  this  date 

# 

ifrcst=2008043018 

# 

#  -  day  when  this  run  is  done  (just  in  output  and  script  being  submitted) 

# 

idtgtod=2 00805011 8 

# 

#  Change  note: 

#  Removed  the  TMPDIR  exclusion. 


###  @  environment  =  COPY  ALL;  ! TMPDIR; 


# 

# 

# 

# 

# 

# 

# 

# 

# 

# 


Script : 
Purpose : 

Author : 

Created: 
Platform: 
Usage : 


ncoda_MERpal . job 

to  setup  and  cycle  the  CODA  analysis  only 

in  the  HYCOM  1/12  deg  Gulf  of  Mexico  grid 

Steve  Lowder,  NRL-CSC 

hacked  by  Jim  Cummings 

1  Mar  2004 

IBM  SP5  at  NAVO 

see  usage ()  function,  below 
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# - 

#  DATA  NOTE: 

#  This  script  copies  data  from  a  mass  storage  server  to 

#  a  local  archive  area  (scratch) .  When  the  analysis 

#  runs  it  operates  in  another  run  directory  (scratch) .  At  the 

#  end  of  the  run,  this  script  does  not  move  the  results  to 

#  a  mass  storage  server.  The  move  to  the  permanent  storage 

#  is  done  with  the  plotting  script  submitted  at  the  end  of  the  script. 

#  - 

# - 

#  FUNCTIONS 

#  Most  of  the  work  in  this  script  is  organized  into 

#  functions  which  are  defined  at  the  beginning  of  the 

#  script.  Search  for  "Begin  script"  to  skip  them. 

#  - 

function  usage  { 

print  ncoda_MERpal . j ob 

print  "usage:  ncoda_MERpal . j ob  start_dtg  num_days  frcst_dtg" 

print  "  start_dtg  is  the  first  day  in  yyyymmddhh  format" 

print  "  num_days  is  the  total  number  of  days  to  cycle" 

print  "  frcst_dtg  is  the  forecast  day  in  yyyymmddhh  format" 

} 


function 

echo 

echo 

echo 

echo 

echo 

echo 

echo 

echo 

echo 

exit 


fatal  { 

I!  II 

II  II 


II 

II 

Script  stopped  on  error  condition  at  "  'date  +%T 

II 


II 


1 


function  setup_paths  { 

# 

#  There  are  three  logical  areas  to  this  script. 

#  1 .  The  mass  storage  areas  where  data  is  permanently  stored. 

#  2.  The  archive  area  on  the  computer  where  mass  storage  data  is 

#  copied  prior  to  a  run. 

#  3.  The  run  area  on  the  computer  where  the  application  does  its 

#  work. 

#  Setup  paths  to  the  mass  storage  server 

# 

#  SAL  change: 

T=0  9 
E=727 

EXPT=expt_7  2 . 7 
REG=MERpalO . 08 
REGl=MERpal 
REG0=GLBa0 . 08 
MSASv=newton 

MSAS_LOCA_MDL=/ scr/ $ { USER} /hycom/ $ { REG} /$ {EXPT } / raw/ nowcast 
MSAS_WORK_TOP=/ u/home/ $ {USER} /hycom/ $ { REGO } / subtopo/ $ { REG1 } 
MSAS_WORK_MDL=/ u/home/ $ {USER} /hycom/ ${ REGO } /means td 
MSAS_WORK_OBS=$ {MSASv} : /u/home/ ooc/data/ncoda/ocnqc 

# 

#  Setup  paths  on  the  computer  in  the  user's  scratch  and  home  directory 
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# 


# 

BIN  DIR 

-  the 

# 

SCRATCH  DIR 

-  the 

# 

MODEL  DIR 

-  the 

# 

DB  DIR 

-  the 

# 

OBS  DIR 

-  the 

# 


location 

location 

location 

location 

location 


of  all  executables 

of  user's  scratch  area 

of  model  specific  static  data 

of  static  database  files  like  dim, land-sea 

of  ocean  data  observations  like  ship,  ssmi 


BIN_DIR=/ u/home/ j  ac/ncoda_mpi/bin 
DDTG=/u/home/ j  ac/ ncoda_dtg/bin/ dtg 

#  BIN_DIR= /u/home/ ooc /mode ls/ncoda/ncoda_mpi /bin 

#  DDTG=/ u/home /ooc /model s/ncoda/ ncoda_dtg/bin/dtg 
SCRATCH_DIR=/scr/$USER 

RUN_DIR=$  SCRATCH_DIR/ $ { REG } / $ { EXPT } 
SAVE_DIR=$SCRATCH_DIR/ $ { REGO } /ncoda/ $ { EXPT } 
RESTART_DIR=$RUN_DIR/ restart 
DB_DIR=$SCRATCH_DIR/ database 
OBS_DIR=$DB_DIR/ ocnqc 


function  setup_dirs  { 

Cd  $SCRATCH_DIR 

export  TMPDIR=$SCRATCH_DIR/tmp 
mkdir  -p  $TMPDIR 


# 

#  Make  the  static  database  dir,  obs,  dir,  and  model  dir  if 

#  they  do  not  exist. 

# 

mkdir  -p  $DB_DIR 
mkdir  -p  $  {  DB_DIR} /dim 
mkdir  -p  $ { DB_DIR} /gdem 


mkdir  -p  $OBS_DIR 

mkdir  -p  $ { OBS_DIR} /beta 

cd  $ { OBS_DIR} /beta 

mkdir  -p  altim 

mkdir  -p  goes 

mkdir  -p  lac 

mkdir  -p  profile 

mkdir  -p  ship 

mkdir  -p  ssmi 

mkdir  -p  ${ OBS_DIR} /gamma 

cd  ${ OBS_DIR} /gamma 

mkdir  -p  altim 

mkdir  -p  goes 

mkdir  -p  lac 

mkdir  -p  profile 

mkdir  -p  ship 

mkdir  -p  ssmi 

mkdir  -p  ${ OBS_DIR} /kappa 

cd  ${ OBS_DIR} /kappa 

mkdir  -p  altim 

mkdir  -p  goes 

mkdir  -p  lac 

mkdir  -p  profile 

mkdir  -p  ship 

mkdir  -p  ssmi 

mkdir  -p  $ { OBS_DIR} /godae 

cd  $ { OBS_DIR} / godae 

mkdir  -p  altim 
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# 

# 

# 

# 

# 

# 


} 


mkdir 

mkdir 

mkdir 

mkdir 


-p  mcsst 
-p  profile 
-p  ship 
-p  ssmi 


Now  remake  the  run  directory  to  be  sure  it  is  clean. 


rm  -rf  $RUN_DIR/analysis 

rm  -rf  $RUN_DIR/ restart 

mkdir  -p  $RUN_DIR 
mkdir  -p  $RUN_DIR/analysis 
mkdir  -p  $RUN_DIR/restart 
mkdir  -p  $RUN_DIR/output 
mkdir  -p  $SAVE_DIR 


function  setup_data  { 

# 

#  Define  local  variables 

# 

typeset  last_dtg 
typeset  -i  n_back  f_frwd 

# 

#  Prepare  the  data  for  the  entire  run.  You  will  need  data  from 

#  three  logical  areas: 

#  1.  Static  data  2.  Observation  data  3.  Model  area  data 

# 

#  Build  or  replenish  the  Model  area  data 

# 

echo  "  " 

echo  "==================================================== 

echo  "  " 


# 

#  Set  the  parameters  for  the  number  of  days  forward 

#  and  backward  to  look  for  obs. 

# 

let  n_back=-10 

let  n_f rwd=2+$ { inumd} 

# 

print  -n  "  Checking  depth  file  (date  time  group  indpendent) " 

Cd  $ { RESTART_DIR} 

#  if  [  !  -r  depths_sfc_000000_000000_lo2525xl841_$ {prev_dtg}_00000000_datafld 
] ;  then 

/bin/cp  $MSAS_WORK_TOP/ depth_$ { REG } _$ { T } .A 
depths_sf c_000000_000000_lo2525xl841_$ {prev_dtg}_00000000__dataf Id 


#  fi 

echo  " .  ready" 

print  -n  "  Checking  initial  model  error  file  (should  be  name  with 

initial  date) " 


Cd  $ { RESTART_DIR} 

#  if  [  !  -r  seahgt_sf c_000000_000000_lo2 52 5x18 41_1 99 9080 900_00000000_modlerr 
] ;  then 

#  rep 

$MSAS_WORK_MDL/ seahgt_s f c_0 00 00 0_00 00 00_lo2  52  5x1 8  4 1_1 9 9  90 8 0  90 0_00 00 00 00_modlerr 

#  rep 

$MSAS_WORK_MDL/ seatmp_s f c_0 00 00 0_00 00 00_lo2  52  5x1 8  4 1_1 9 9  90 8 0  90 0_00 00 00 00_modlerr 

#  rep 

$MSAS_WORK_MDL/ seatmp_pre_0 00 00 0_002 5 00_lo2  52  5x1 8  4 1_1 9 9  90 8 0  90 0_00 00 00 00_modlerr 

#  rep 

$MSAS_WORK_MDL/ salint_pre_0 00 00 0_002 5 00_lo2  52  5x1 8  4 1_1 9 9  90 8 0  90 0_00 00 00 00_modlerr 
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#  rep 

$MSAS_WORK_MDL/uucurr_pre_000000_002500_lo2  52  5x18  41_1 99  9080  900_00000000_modlerr 

#  rep 

$MSAS_WORK_MDL/ vveurr_pre_0 00 00 0_002 5 00_lo2  52  5x1 8  4 1_1 9 9  90 8 0  90 0_00 00 00 00_modlerr 

#  fi 

#  echo  " .  ready" 

FCST=${ifrcst} 

echo  'FORECAST  DATE  ’${FCST} 


print  -n  "  Checking  first  guess  fields  valid  tau  024" 

Cd  $ { RESTART_DIR} 

#  if  [  !  -r  seatmp_sf c_000000_000000_lo2  52  5x18  41_$ { FCST }_002  4  0000_f cstf Id  ]; 
then 

#  rep  $MSAS_WORK_MDL/hycom2dt_000000_$ { FCST }_0024 0000 .A 
seatmp_sf c_000000_000000_lo2525xl841_$ { FCST }_0024 0000_f cstf Id 

/ bi n / cp  $MSAS_LOCA_MDL/hycom2 dt_0 0  0  0  0  0_$ { FC  ST } _0  0  2  4  0  0  0  0 . A 
seatmp^sf c_000000_000000_lo2525xl841_$ { FCST }_0024 0000_f cstf Id 

#  fi 

#  if  [  !  -r  seahgt_sf c_000000_000000_lo2  52  5x18  41_$ { FCST }_002  4  0000_f cstf Id  ]; 
then 

#  rep  $MSAS_WORK_MDL/hycom2dSf sd_000000_$ { FCST }_0 02 4  0000 . A 
seahgt_sfc_000000_000000_lo2525xl841_$ {FCST }_00240000_fcstfld 

/bin/cp  $MSAS_LOCA_MDL/hycom2dSf sd_000000_$ {FCST }_0 02 40000 . A 
seahgt_sf c_000000_000000_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 
#for  tau  in  00030000  00060000  00090000  00120000  00150000  00180000  00210000 
00240000 
#do 

#  /bin/cp  $MSAS_LOCA_MDL/hycom2dSf sd_000000_$ { FCST }_$ { tau } .A 
seahgt_sf c_000000_000000_lo2525xl841_$ { FCST }_$ { tau}_f cstf Id 

#done 

#  fi 

#  if  [  !  -r  seaice_sf c_000000_000000_lo252 5x18 41_$ { FCST }_002 4  0000_f cstf Id  ]; 

then 

#  rep  $MSAS_WORK_MDL/hycom2dice_000000_$ { FCST }_0 02 4  0000 . A 
seaice_sfc_000000_000000_lo2525xl841_$ {FCST }_00240000_fcstfld 

/bin/cp  $MSAS_LOCA_MDL/hycom2 di ce_0 00000_${FCST}_00240000.A 
seaice_sf c_000000_000000_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 

#  fi 

#  if  [  !  -r  mixlyr_sf c_000000_000000_lo2 52 5x18 41_$ { FCST }_002 4  0000_f cstf Id  ]; 

then 

#  rep  $MSAS_WORK_MDL/hycom2dmix_000000_$ {FCST }_0 02 40000 . A 
mixlyr_sf c_000000_000000_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 

/bin/cp  $MSAS_LOCA_MDL/hycom2dmix_000000_$ {FCST }_0024 0000 .A 
mixlyr_sf c_000000_000000_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 

#  fi 

#  if  [  !  -r  seahgt_sf c_000000_000000_lo2 52 5x18 41_$ { FCST }_00000000_meanf Id  ]; 

then 

/bin/cp  $MSAS_WORK_MDL/RS_0 94_$ { T }_$ { REG1 } 0 . 08 . A 
seahgt_sf c_000000_000000_lo2525xl841_$ { FCST }_00 0000 00_meanf Id 

#  rep  $MSAS_WORK_MDL/ 056_archMNA. 0009_0013_fsd. A 

seahgt_sf c_000000_000000_lo2525xl841_$ { FCST }_00 0000 00_meanf Id 

#  rep  $MSAS_WORK_MDL/micomecmwf . year04_05_$ { REG } . A 

seahgt_sf c_000000_000000_lo2525xl841_$ { FCST }_00 0000 00_meanf Id 

#  rep  $MSAS_WORK_MDL/ zero_0 08_hycom . A 

seahgt^sf c_000000_000000_lo2525xl841_$ { FCST }_00 0000 00_meanf Id 

#  fi 

#  if  [  !  -r  densty_pre_000000_002500_lo2 52 5x18 41_$ { FCST }_002 4  0000_f cstf Id  ]; 

then 

##  rep  $MSAS_WORK_MDL/hycom3dden_002500_$ {FCST }_0 02 40000 . A 

densty_pre_000000_002  500_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 

#  /bin/cp  $MSAS_LOCA_MDL/hycom3dden_002500_$ { FCST }_0 02 4  0000 . A 
densty_pre_000000_002  500_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 
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#  fi 

#  if  [  !  -r  seatmp_pre_000000_002500_lo2  52  5x18  41_$ { FCST }_002  4  0000_f cstf Id  ]; 
then 

#  rep  $MSAS_WORK_MDL/hycom3dt_002500_$ {FCST }_0024  0000 .A 

seatmp_pre_000000_002500_lo2525xl841_$ { FCST }_0024 0000_f cstf Id 

/bin/cp  $MSAS_LOCA_MDL/hycom3dt_002  500_$ {FCST }_0024 0000 .A 
seatmp_pre_000000_002500_lo2525xl841_$ { FCST }_00240000_fcstfld 

#  fi 

#  if  [  !  -r  salint_pre_000000_002500_lo2 52 5x18 41_$ { FCST }_002 4  0000_f cstf Id  ]; 

then 

#  rep  $MSAS_WORK_MDL/hycom3ds_002500_$ { FCST }_00240000 .A 

salint_pre_000000_002  500_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 

/bin/cp  $MSAS_LOCA_MDL/hycom3ds_002 500_$ {FCST }_0024  0000 .A 
salint_pre_000000_002  500_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 

#  fi 

#  if  [  !  -r  uucurr_pre_000000_002500_lo2  52  5x18  41_$ { FCST }_002 4  0000_f cstf Id  ]; 

then 

#  rep  $MSAS_WORK_MDL/hycom3du_002  500_$ {FCST }_0024 0000 .A 

uucurr_pre_000000_002500_lo2525xl841_$ { FCST }_0024 0000_f cstf Id 

/bin/cp  $MSAS_LOCA_MDL/hycom3du_002 500_$ {FCST }_0024 0000 .A 
uucurr_pre_000000_002500_lo2525xl841_$ { FCST }_0024 0000_f cstf Id 

#  fi 

#  if  [  !  -r  vvcurr_pre_000000_002500_lo2 52 5x18 41_$ { FCST }_002 4  0000_f cstf Id  ]; 

then 

#  rep  $MSAS_WORK_MDL/hycom3dv_002500_${FCST}_00240000.A 

vvcurr_pre_000000_002  500_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 

/bin/cp  $MSAS_LOCA_MDL/hycom3dv_002 500_$ {FCST }_0024 0000 .A 
vvcurr_pre_000000_002  500_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 

#  fi 

if  [  !  -r  seahgt_sf c_000000_000000_lo2525xl 84 1_$ { FCST }_0024 0000_f cstf Id  ]; 

then 

echo  "seahgt  sfc  forecast  does  not  exist  for  "${FCST} 

exit 

fi 

if  [  !  -r  seatmp_sf c_000000_000000_lo2525xl 84 1_$ { FCST }_0024 0000_f cstf Id  ]; 

then 

echo  "seatmp  sfc  forecast  does  not  exist  for  "${FCST} 

exit 

fi 

if  [  !  -r  seatmp_pre_000000_002 500_lo2525xl 84 1_$ { FCST }_0024 0000_f cstf Id  ]; 

then 

echo  "seatmp  forecast  does  not  exist  for  "${FCST} 

exit 

fi 

if  [  !  -r  salint_pre_000000_002 500_lo2525xl 84 1_$ {FCST }_0024 0000_f cstf Id  ]; 

then 

echo  "salint  forecast  does  not  exist  for  "${FCST} 

exit 

fi 

if  [  !  -r  uucurr_pre_000000_002 500_lo2525xl 84 1_$ { FCST }_0024 0000_f cstf Id  ]; 

then 

echo  "uucurr  forecast  does  not  exist  for  "${FCST} 

exit 

fi 

if  [  !  -r  vvcurr_pre_000000_002 500_lo2525xl 84 1_$ { FCST }_0024 0000_f cstf Id  ]; 

then 

echo  "vveurr  forecast  does  not  exist  for  "${FCST} 


exit 

fi 

echo  " .  ready" 

print  -n  "  Checking  binary  files  " 
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if  [  !  -r  $BIN_DIR/ncoda_prep  ] ;  then 
echo  "Error:  ncoda_prep  not  found" 
fi 

if  [  !  -r  $BIN_DIR/ncoda  ] ;  then 
echo  "Error:  ncoda  not  found" 
fi 

if  [  !  -r  $BIN_DIR/ncoda_post  ] ;  then 
echo  "Error:  ncoda_post  not  found" 
fi 

echo  " .  ready" 

echo  "  " 
echo  "  " 

echo  "================================: 

echo  "  " 

echo  "  Data  preparation  complete:" 
echo  "  " 

echo  "================================: 


function  setup_processors  { 

# 

#  Set  OpenMP  environment 

#  Set  number  of  OpenMP  processors 

# 

#  Note:  The  pre  and  post  analysis  use  OpenMP  and  the 

#  analysis  uses  MPI  only. 

# 

#  Set  number  of  MPI  processors 

#  MPICMD= 'mpirun . lsf  ' 

#  MPICMD='poe  ' 
nprc=4 1 6 

#  Set  number  of  OpenMP  processors 

mprc=l 

# 

#  Set  OpenMP  environment 

# 

( (SLAVE=512*1024*1024 ) ) 
export  OMP_NUM_THREADS=$ {mprc } 

#  IBM 

export  XLSMPOPTS="stack=$SLAVE" 

#  XT5 

#  export  MPSTKZ=512M 

export  OMP_SCHEDULE="DYNAMIC, 1 " 

# 

#  export  OMPCMD= ' timex  -p  -mt  ' 

#  export  MPICMD= ' timex  -p  -mt  'poe 

#  export  MPICMD= ' timex  -p  -mt  'mpirun. lsf 

#  XT5 

export  OMPCMD= ' aprun  -n  1  -d  ' $ { OMP_NUM_THREADS } 
export  MPICMD= ' aprun  -n  '${nprc} 


echo  "  " 

echo  "##  PE  Environment  Overview  {" 
printenv  |  grep  -e  ,AOMP_'  -e  ' AMP '  -e 
'  AMPI  ' 


echo  "##  }" 
echo  "" 

# 

#  Print  the  hard  and  soft  limits. 

# 


DSM '  -e  ' A CHUNK '  - 


'  A  PAGE '  -e 
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ulimit  -aH 
ulimit  -aS 


function  setup_namelists  { 

# 

#  Setup  the  namelist  files 

# 

rm  -f  odsetnl 
rm  -f  gridnl 
rm  -f  oanl 

# 

#  data  path  settings 

# 

.  /u/home/$ {USER} /$ {REGO } /ncoda/$ { EXPT } /$ { REG1 }_$ {E } . odsetnl 

# 

#  grid  definition  namelist  settings 

#  kko  is  the  number  of  vertical  levels;  levels  are  defined  in  oanl 

# 

.  / u/home/ $ {USER} /$ { REGO } /ncoda/$ { EXPT } /$ { REG1 }_$ { E } .gridnl 

# 

#  ocean  analysis  namelist  settings 

#  for  CH  assimilation  of  altimeter  SSH,  set  direct  =  .true. 

#  for  MODAS  assimilation  of  altimeter  SSH,  set  modas  =  .true. 

#  for  SST  assimilation,  set  st_asm  =  .true. 

#  for  potential  temperature  analysis,  set  pt_anl  =  .true. 

#  (assumes  potential  temperature  background  fields) 

# 

.  / u/home/ $ {USER} /$ { REGO } /ncoda/$ { EXPT } /$ { REG1 }_$ { E } .oanl 

# 

echo  "  " 

} 


function  run_analysis  { 
typeset  analysis_dtg 
typeset  -i  cycle^index 

analysis_dtg=$ { idtg} 
cycle_index=$ { inumd} 

# 

#  Remove  the  work  files  created  by  the  analysis 

# 

/bin/rm  -f  pout* 


echo  "  Start  Ocean  Prep  2D  at  "  'date  +%T'  »  poutl 
echo  "  "  »  poutl 

${OMPCMD}  $BIN_DIR/ncoda  prep  2D  gridnl  $analysis  dtg  >>  poutl 


if  [  $?  -ne 
echo 
echo 
echo 
echo 
echo 
echo 
fatal 
else 
echo 
echo 
echo 
echo 
echo 
echo 
fi 


0  ] ;  then 
»  poutl 
»  poutl 

Error;  Ocean  Prep  2D  failed  at 
Check  poutl"  »  poutl 
»  poutl 
»  poutl 


»  poutl 
»  poutl 
End  Ocean 
Completed 
»  poutl 
»  poutl 


Prep  2D  at  " 
successfully" 


'date  +%T '  »  poutl 


date  +%T ' 
»  poutl 


>>  poutl 
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echo  "  Start  Ocean  Analysis  2D  at  "  'date  +%T'  »  pout2 
echo  "  "  >>  pout2 

${MPICMD}  $BIN_DIR/ncoda  2D  gridnl  $analysis_dtg  »  pout2 
if  [  $?  -ne  0  ] ;  then 
echo  "  "  »  pout2 
echo  "  "  »  pout2 

echo  "  Error:  Ocean  Analysis  2D  failed  at  "  'date  +%T'  »  pout2 
echo  "  Check  pout2"  »  pout2 
echo  "  "  »  pout2 
echo  "  "  »  pout2 
fatal 
else 

echo  "  "  »  pout2 
echo  "  "  »  pout2 

echo  "  End  Ocean  Analysis  2D  at  "  'date  +%T'  >>  pout2 

echo  "  Completed  successfully"  »  pout2 

echo  "  "  »  pout2 
echo  "  "  »  pout2 


echo  "  Start  Ocean  Post  2D  at  "  'date  +%T'  »  pout3 
echo  "  "  >>  pout3 

${OMPCMD}  $BIN_DIR/ncoda_post  2D  gridnl  $analysis_dtg  >>  pout3 
if  [  $?  -ne  0  ] ;  then 
echo  "  "  »  pout3 
echo  "  "  »  pout3 

echo  "  Error:  Ocean  Post  2D  failed  at  "  'date  +%T'  >>  pout3 
echo  "  Check  pout3"  »  pout3 
echo  "  "  »  pout3 
echo  "  "  »  pout3 
fatal 
else 

echo  "  "  »  pout3 
echo  "  "  »  pout3 

echo  "  End  Ocean  Post  2D  at  "  'date  +%T'  >>  pout3 
echo  "  Completed  successfully"  >>  pout3 
echo  "  "  »  pout3 
echo  "  "  »  pout3 


echo  "  Start  Ocean  Prep  3D  at  "  'date  +%T'  »  pout4 
echo  "  "  »  pout4 

${OMPCMD}  $BIN_DIR/ncoda_prep  3D  gridnl  $analysis_dtg  >>  pout4 
if  [  $?  -ne  0  ] ;  then 
echo  "  "  »  pout4 
echo  "  "  »  pout4 

echo  "  Error:  Ocean  Prep  3D  failed  at  "  'date  +%T'  >>  pout4 

echo  "  Check  pout4"  »  pout4 

echo  "  "  »  pout4 
echo  "  "  »  pout4 
fatal 
else 

echo  "  "  »  pout4 
echo  "  "  »  pout4 

echo  "  End  Ocean  Prep  3D  at  "  'date  +%T'  >>  pout4 
echo  "  Completed  successfully"  >>  pout4 
echo  "  "  »  pout4 
echo  "  "  »  pout4 


echo  "  Start  Ocean  Analysis  3D  at  "  'date  +%T'  »  pout5 
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echo  "  "  >>  pout5 

${MPICMD}  $BIN_DIR/ncoda  3D  gridnl  $analysis_dtg  »  pout5 
if  [  $?  -ne  0  ] ;  then 

if  ( (  $cycle_index  ==  1  ) ) ;  then 
echo  "  "  »  pout5 
echo  "  "  »  pout5 

echo  "  Error:  Ocean  Analysis  3D  failed  at  "  'date  +%T'  >>  pout5 

echo  "  Check  pout5"  »  pout5 

echo  "  This  is  not  a  fatal  condition,  continuing  ..."  »  pout5 
echo  "  "  »  pout5 
echo  "  " 
else 

echo  "  "  »  pout5 
echo  "  "  »  pout5 

echo  "  Error:  Ocean  Analysis  3D  failed  at  "  'date  +%T'  >>  pout5 

echo  "  Check  pout5"  »  pout5 

echo  "  "  »  pout5 
echo  "  "  »  pout5 
fatal 
fi 

else 

echo  "  "  »  pout5 
echo  "  "  »  pout5 

echo  "  End  Ocean  Analysis  3D  at  "  'date  +%T'  >>  pout5 

echo  "  Completed  successfully"  »  pout5 

echo  "  "  »  pout5 
echo  "  "  »  pout5 


echo  "  Start  Ocean  Post  3D  at  "  'date  +%T'  »  pout6 
echo  "  "  >>  pout 6 

${OMPCMD}  $BIN_DIR/ncoda_post  3D  gridnl  $analysis_dtg  >>  pout6 
if  [  $?  -ne  0  ] ;  then 

if  ( (  $cycle_index  ==  1  ) ) ;  then 
echo  "  "  »  pout6 
echo  "  "  »  pout6 

echo  "  Error:  Ocean  Post  3D  failed  at  "  'date  +%T'  >>  pout6 

echo  "  Check  pout6"  »  pout6 

echo  "  This  is  not  a  fatal  condition,  continuing  ..."  »  pout6 
echo  "  "  »  pout6 
echo  "  "  »  pout6 
else 

echo  "  "  »  pout6 
echo  "  "  »  pout6 

echo  "  Error:  Ocean  Post  3D  failed  at  "  'date  +%T'  >>  pout6 

echo  "  Check  pout6"  »  pout6 

echo  "  "  »  pout6 
echo  "  "  »  pout6 
fatal 
fi 

else 

echo  "  "  »  pout6 
echo  "  "  »  pout6 

echo  "  End  Ocean  Post  3D  at  "  'date  +%T'  >>  pout6 
echo  "  Completed  successfully"  »  pout6 
echo  ""  >>  pout 6 
echo  ""  >>  pout6 
fi 

# 

#  rename  diagnostic  files  (creation  controlled  by  debug  options  in  oanl) 

# 

mv  fort . 25  $ { SAVE_DIR} /$ { REG1 } _ $ { E } _ $ { idtgtod}_$ {analysis_dtg} . rpr 

mv  fort . 27  $ { SAVE_DIR} / $ { REG1 } _ $ { E } _ $ { idtgtod}_$ {analysis_dtg} . vfy 
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# 

# 

# 


mv 

fort . 

31 

mv 

fort . 

32 

mv 

fort . 

33 

mv 

fort . 

34 

mv 

fort . 

35 

mv 

fort . 

36 

mv 

fort . 

37 

mv 

fort . 

38 

mv 

fort . 

39 

$ { SAVE_DIR} /$ { REG1 } _ $ { E } _ $ { idtgtod }_$ {analysis_dtg} . syn 

$ { SAVE_DIR} / $ { REG1 } _ $ { E }_$ { idtgtod} _$ { analysi s_dtg } . rej 

$ { SAVE_DIR} /$ {REG1 } _ $ { E } _ $ { idtgtod }_$ { analysi s_dtg} .prf 

$ {SAVE_DIR} /$ { REG1 } _$ { E }_$ { idtgtod} _$ { analysi s_dtg } .gpt 

${SAVE_DIR} / $ { REG1 } _ $ { E }_$ { idtgtod }_$ { analysi s_dtg} .err 

$ { SAVE_DIR} /$ { REG1 } _ $ { E } _ $ { idtgtod }_$ { analysi s_dtg} .mvo 

$ { SAVE_DIR} /$ {REG1 } _ $ { E } _ $ { idtgtod }_$ { analysi s_dtg} . drc 

${SAVE_DIR} /${REG1 }_$ { E }  $ { idtgtod} _$ {analysi s_dtg} . lyp 
${SAVE_DIR} / $ { REG1 }_$ { E }_$ { idtgtod }_$ { analysi s_dtg} . fix 


combine  program  outputs  into  single  file  and  finish 


/bin/rm  $analysis_dtg . out 

cat  poutl  pout2  pout3  pout4  pout5  pout6  >>  $analysis_dtg . out 

echo  " -  odsetnl  - "  »  $analysis_dtg . out 

cat  odsetnl  »  $analysis_dtg. out 

echo  " -  gridnl  - "  >>  $analysis_dtg. out 

cat  gridnl  »  $analysis_dtg . out 

echo  " - oanl - "  >>  $analysis_dtg.  out 

cat  oanl  »  $analysis_dtg . out 

echo  "  Analysis  ended  at  "  'date  +%T' 

echo  "  Output  file:  $analysis_dtg . out  in  "  'pwd' 

echo  "  "  'date  -u  +'%Y%m%d'00' 

echo  "  " 

echo  "  " 

echo  "====================================================" 

mv  $analysis_dtg. out  $ { SAVE_DIR} / $ { REG1 } $ { E } $ { idtgtod}_$ { analysis_dtg } . out 


function  cleanup  { 

cd  $RUN_DIR/analysis 

#  /bin/rm  -f  coda.*obs& 

#  /bin/ls  |  egrep  ' AMVOI '  |  xargs  /bin/rm 

/bin/rm  MVOI*& 

#  rm  -f  pout* 

} 


# 

# - 

# - 

# 

#  Begin  script  execution,  check  command  line  arguments 

# 

# - 

# - 

# 

#if  [  "$#"  -ne  3  ];  then 

#  clear 

#  usage 

#  print  ERROR:  incorrect  number  of  input  args 

#  exit  1 
#else 

FIRST_DTG="$ { idtg} " 

NUM_DAYS=" $ { inumd} " 

FCST="${ifrcst}" 

#fi 

if  [[  $ { FIRST_DTG }  !=  +([0-9])  ]];  then 
clear 
usage 

print  ERROR:  arg  1  is  not  an  integer 
exit  1 
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fi 

if  [[  ${NUM_DAYS]  !=  +([0-9])  ]];  then 
clear 
usage 

print  ERROR:  arg  2  is  not  an  integer 
exit  1 
fi 

if  [[  $ { FCST }  ! =  +([0-9])  ]];  then 
clear 
usage 

print  ERROR:  arg  3  is  not  an  integer 
exit  1 
fi 

if  [  $ { #FIRST_DTG }  -ne  10  ];  then 
clear 
usage 

print  ERROR:  incorrect  format  for  arg  1 
exit  1 
fi 

if  [  $ { #FCST }  -ne  10  ];  then 
clear 
usage 

print  ERROR:  incorrect  format  for  arg  3 
exit  1 
fi 

# 

#  The  batch  process  may  set  file  creation  mask  to  600,  change 

#  this  with  the  mask. 

# 

umask  022 

# 

#  Define  some  useful  variables  for  the  model  name,  dimensions, 

#  and  number  of  hours  in  update  cycle. 

# 

MODEL=hycom 
UPDATE_CYCLE=24 
setup_paths  $ [MODEL] 

# 

#  Setup  date-time-groups  for  the  script 

# 

let  "day_index=$ { NUM_DAYS }  -  1" 
export  CRDATE=$ [FIRST_DTG: -19 99060600} 

prev_dtg= ' $DDTG  -d  -1  2>  /dev/null' 

#  prev_dtg= ' $DDTG  -d  -7  2>  /dev/null' 

LAS T_DTG= ' $ DDTG  -d  ${day_index}  2>  /dev/null' 
curr_dtg= ' $DDTG  2>  /dev/null' 


echo 

=  =  =  = 

— 

— 

echo 

II 

II 

echo 

II 

II 

echo 

"CODA 

analysis  processor (ncoda 

$ [REG1 } .job) " 

echo 

II 

The  current  date/time  is" 

$ (  date) 

echo 

II 

II 

echo 

II 

II 

echo 

II 

Dates :  " 

echo 

II 

First  cycle  date 

:  $FIRST  DTG" 

echo 

II 

Last  cycle  date 

:  $LAST  DTG" 
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echo  "  Forecast  date  :  $FCST" 

echo  "  " 
echo  "  " 

echo  "  Path  setting:" 

echo  "  Binary  locations:  $BIN_DIR" 

echo  "  Scratch  directory:  $SCRATCH_DIR" 

#  echo  "  Model  directory:  $MODEL_DIR" 

echo  "  Restart  directory:  $RESTART_DIR" 

echo  "  CODA  database  directory:  $DB_DIR  " 

echo  "  CODA  obs  data  directory:  $OBS_DIR  " 

echo  "  " 
echo  "  " 

echo  "====================================================" 


setup_dirs 

# 

#  Clean  the  temporary  directory  structure 

# 

if  [[  !  -e  $RUN_DIR  ]];  then 

echo  "The  run  directory  is  missing,  can  not  proceed:  $RUN_DIR" 
exit  1 
fi 

if  [ [  !  -e  $RUN_DIR/restart  ] ] ;  then 

echo  "The  Restart  directory  is  missing,  can  not  proceed:  $RUN_DIR/ restart" 
exit  1 
fi 

cd  $RUN_DIR 

setup_data  $ { FIRST_DTG}  ${NUM_DAYS}  ${FCST} 
setup_processors 

# 

#  Execute  the  coda  mvoi  programs 

# 

echo  "  " 
echo  "  " 

echo  "Start  Analysis" 

echo  "  Directories  created" 

echo  "  Date-time-group :  $curr_dtg" 

echo  "  MPI  processes:  $nprc" 

echo  "  OpenMP  threads:  $mprc" 

echo  "  Initial  directory  "  'pwd' 

echo  "  " 

cd  $RUN_DIR/analysis 
setup_namelists 

# 

#  loop  over  update  cycles 

# 

i=l 

h=0 

num_cycles=$ { NUM_DAYS } 

#while  ( (  $i  <=  $num_cycles  ) )  do 

# 

#  set  analysis  dtg 

# 

curr_dtg= ' $DDTG  -h  $h  2>  /dev/null' 
echo  'UPDATE  CYCLE  :  ' $UPDATE  CYCLE 
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echo  'curr_dtg  :  '${curr_dtg}  $ { i } 
idtg=${idtg} 

#  execute  ncoda  mpi  mvoi  programs 

run_analysis  ${idtg}  ${i} 

#  run_analysis  ${curr_dtg}  $ { i } 

#  update  counters 

#  (  (  i  =  i  +  1  )  ) 

#  (  (  h  =  h  +  $UPDATE_CYCLE  )  ) 

#done 

♦cleanup 

# 

cd  /scr/$ {USER} /hycom/ $ { REGO } /$ { EXPT } /logs 
touch  done_$ { REG1 }_$ { idtgtod}_$ { curr_dtg} 

/u/home/wallcraf /bin/q_navo  $ { E } ncoda2arch_$ { idtgtod}_$ { curr_dtg} . csh 

# 

/u/home/wallcraf /bin/ q_navo  $ {E }plot_$ { REG1 }_$ { idtgtod}_$ { curr_dtg} . csh 
echo  "" 

echo  "====================================================" 

exit  0 
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APPENDIX  E 

Script  727ncoda2arch_${idtgtod}18_${idtg}18.com 

This  script  converts  the  NCODA  analysis  on  z-levels  to  HYCOM  vertical  coordinates. 


# ! /bin/ksh 
#PBS  -N  999pbs 
#PBS  -j  oe 

#PBS  -o  999pbsXX.log 
#PBS  -1  mppwidth=416 
#PBS  -1  mppnppn=8 
#PBS  -1  walltime=4 : 00 : 00 
#PBS  -W  umask=027 
#PBS  -A  NRLSS018 
#PBS  -q  internal3d 
# 

#BSUB  -J  727ncoda_MERpal . job 
#BSUB  -n  192 
#BSUB  -a  poe 
#BSUB  -W  3:30 

#BSUB  -R  "rusage [ntbl_windows=32 ] span [ptile=l 6 ] " 

#BSUB  -q  challenge 
#BSUB  -P  NRLSSC3 J 
#####BSUB  -q  internal3d 
#####BSUB  -P  NAVOSHYC 
#BSUB  -e  err_MERpal . % J 
# 

#  -  set  days  to  run 

# 

inumd=l 

# 

#  -  analysis  date 

# 

idtg=2008050118 

# 

#  -  forecast  from  this  date 

# 

ifrcst=2008043018 

# 

#  -  day  when  this  run  is  done  (just  in  output  and  script  being  submitted) 

# 

idtgtod=2 00805011 8 

# 

#  Change  note: 

#  Removed  the  TMPDIR  exclusion. 


###  @  environment  =  COPY  ALL;  ! TMPDIR; 


# 

# 

# 

# 

# 

# 

# 

# 

# 

# 


Script : 
Purpose : 

Author : 

Created: 
Platform: 
Usage : 


ncoda^MERpal . job 

to  setup  and  cycle  the  CODA  analysis  only 

in  the  HYCOM  1/12  deg  Gulf  of  Mexico  grid 

Steve  Lowder,  NRL-CSC 

hacked  by  Jim  Cummings 

1  Mar  2004 

IBM  SP5  at  NAVO 

see  usage ()  function,  below 
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# - 

#  DATA  NOTE: 

#  This  script  copies  data  from  a  mass  storage  server  to 

#  a  local  archive  area  (scratch) .  When  the  analysis 

#  runs  it  operates  in  another  run  directory  (scratch) .  At  the 

#  end  of  the  run,  this  script  does  not  move  the  results  to 

#  a  mass  storage  server.  The  move  to  the  permanent  storage 

#  is  done  with  the  plotting  script  submitted  at  the  end  of  the  script. 

#  - 

# - 

#  FUNCTIONS 

#  Most  of  the  work  in  this  script  is  organized  into 

#  functions  which  are  defined  at  the  beginning  of  the 

#  script.  Search  for  "Begin  script"  to  skip  them. 

#  - 

function  usage  { 

print  ncoda_MERpal . j ob 

print  "usage:  ncoda_MERpal . j ob  start_dtg  num_days  frcst_dtg" 

print  "  start_dtg  is  the  first  day  in  yyyymmddhh  format" 

print  "  num_days  is  the  total  number  of  days  to  cycle" 

print  "  frcst_dtg  is  the  forecast  day  in  yyyymmddhh  format" 

} 


function 

echo 

echo 

echo 

echo 

echo 

echo 

echo 

echo 

echo 

exit 


fatal  { 

II  II 

II  II 


II 

II 

Script  stopped  on  error  condition  at  "  'date  +%T 

II 


II 


1 


function  setup_paths  { 

# 

#  There  are  three  logical  areas  to  this  script. 

#  1 .  The  mass  storage  areas  where  data  is  permanently  stored. 

#  2.  The  archive  area  on  the  computer  where  mass  storage  data  is 

#  copied  prior  to  a  run. 

#  3.  The  run  area  on  the  computer  where  the  application  does  its 

#  work. 

#  Setup  paths  to  the  mass  storage  server 

# 

#  SAL  change: 

T=0  9 
E=727 

EXPT=expt_7  2 . 7 
REG=MERpalO . 08 
REGl=MERpal 
REG0=GLBa0 . 08 
MSASv=newton 

MSAS_LOCA_MDL=/scr/$ { USER} /hycom/${ REG} /${EXPT} /raw/nowcast 
MSAS_WORK_TOP=/ u/home/ $ {USER} /hycom/ $ { REGO } / subtopo/ $ { REG1 } 
MSAS_WORK_MDL=/ u/home/ $ {USER} /hycom/$ {REGO } /means td 
MSAS_WORK_OBS=$ {MSASv} : /u/home/ ooc/data/ncoda/ocnqc 

# 
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#  Setup  paths 

# 

#  BIN_DIR 

#  SCRATCH_DIR 

#  M0DEL_DIR 

#  DB_DIR 

#  0BS_DIR 

# 


on  the  computer  in  the  user's  scratch  and  home  directory 

-  the  location  of  all  executables 

-  the  location  of  user's  scratch  area 

-  the  location  of  model  specific  static  data 

-  the  location  of  static  database  files  like  dim, land-sea 

-  the  location  of  ocean  data  observations  like  ship,  ssmi 


BIN_DIR=/ u/ home / j  ac/ncoda_mpi /bin 
DDTG=/u/home/ j ac/ncoda_dtg/bin/ dtg 

#  BIN_DIR=/u/home/ ooc/models/ncoda/ncoda_mpi/bin 

#  DDTG=/ u/home/ooc /model s/ncoda/ ncoda_dtg/bin/ dtg 
SCRATCH_DIR=/scr/$USER 
RUN_DIR=$SCRATCH_DIR/ $ { REG } / $ { EXPT } 
SAVE_DIR=$SCRATCH_DIR/ $ { REGO } /ncoda/$ {EXPT } 
RESTART_DIR=$RUN_DIR/ restart 
DB_DIR=$SCRATCH_DIR/ database 
OBS_DIR=$DB_DIR/ ocnqc 


function  setup_dirs  { 

Cd  $SCRATCH_DIR 

export  TMPDIR=$SCRATCH_DIR/tmp 
mkdir  -p  $TMPDIR 


# 

#  Make  the  static  database  dir,  obs,  dir,  and  model  dir  if 

#  they  do  not  exist. 

# 

mkdir  -p  $DB_DIR 
mkdir  -p  $  {  DB_DIR} /dim 
mkdir  -p  $ { DB_DIR} /gdem 


mkdir  -p  $OBS_DIR 

mkdir  -p  $ { OBS_DIR} /beta 

cd  ${OBS_DIR} /beta 

mkdir  -p  altim 

mkdir  -p  goes 

mkdir  -p  lac 

mkdir  -p  profile 

mkdir  -p  ship 

mkdir  -p  ssmi 

mkdir  -p  ${ OBS_DIR} /gamma 

cd  ${ OBS_DIR} /gamma 

mkdir  -p  altim 

mkdir  -p  goes 

mkdir  -p  lac 

mkdir  -p  profile 

mkdir  -p  ship 

mkdir  -p  ssmi 

mkdir  -p  ${ OBS_DIR} /kappa 

cd  ${OBS_DIR} /kappa 

mkdir  -p  altim 

mkdir  -p  goes 

mkdir  -p  lac 

mkdir  -p  profile 

mkdir  -p  ship 

mkdir  -p  ssmi 

mkdir  -p  $ { OBS_DIR} /godae 

cd  $ { OBS_DIR} / godae 
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# 

# 

# 

# 

# 

# 


} 


mkdir 

mkdir 

mkdir 

mkdir 

mkdir 


-p  altim 
-p  mcsst 
-p  profile 
-p  ship 
-p  ssmi 


Now  remake  the  run  directory  to  be  sure  it  is  clean. 


rm  -rf  $RUN_DIR/analysis 

rm  -rf  $RUN_DIR/ restart 

mkdir  -p  $RUN_DIR 
mkdir  -p  $RUN_DIR/analysis 
mkdir  -p  $RUN_DIR/restart 
mkdir  -p  $RUN_DIR/output 
mkdir  -p  $SAVE_DIR 


function  setup^data  { 

# 

#  Define  local  variables 

# 

typeset  last_dtg 
typeset  -i  n_back  f_frwd 

# 

#  Prepare  the  data  for  the  entire  run.  You  will  need  data  from 

#  three  logical  areas: 

#  1.  Static  data  2.  Observation  data  3.  Model  area  data 

# 

#  Build  or  replenish  the  Model  area  data 

# 

echo  "  " 

echo  "==================================================== 

echo  "  " 


# 

#  Set  the  parameters  for  the  number  of  days  forward 

#  and  backward  to  look  for  obs. 

# 

let  n_back=-10 

let  n_f rwd=2+$ { inumd} 

# 

print  -n  "  Checking  depth  file  (date  time  group  indpendent) " 

Cd  $ { RESTART_DIR} 

#  if  [  !  -r  depths_sfc_000000_000000_lo2525xl841_$ {prev_dtg}_00000000_datafld 
] ;  then 

/bin/cp  $MSAS_WORK_TOP/ depth_$ { REG } _$ { T } .A 
depths_sf c_000000_000000_lo2525xl841_$ {prev_dtg}_00000000_dataf Id 


#  fi 

echo  " .  ready" 

print  -n  "  Checking  initial  model  error  file  (should  be  name  with 

initial  date) " 


Cd  $ { RESTART_DIR} 

#  if  [  !  -r  seahgt_sf c_000000_000000_lo2 52 5x18 41_1 99 9080 900_00000000_modlerr 
] ;  then 

#  rep 

$MSAS_WORK_MDL/ seahgt_s f c_0 00 00 0_00 00 00_1 o2  52  5x1 8  4 1_1 9 9  90 8 0  90 0_00 00 00 00_modlerr 

#  rep 

$MSAS_WORK_MDL/ seatmp_s f c_0 00 00 0_00 00 00_lo2  52  5x1 8  4 1_1 9 9  90 8 0  90 0_00 00 00 00_modlerr 

#  rep 

$MSAS_WORK_MDL/ seatmp_pre_0 00 00 0_002 5 00_lo2  52  5x1 8  4 1_1 9 9  90 8 0  90 0_00 00 00 00_modlerr 

#  rep 

$MSAS_WORK_MDL/ salint_pre_0 00 00 0_002 5 00_lo2  52  5x1 8  4 1_1 9 9  90 80  90 0_00 00 00 00_modlerr 
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#  rep 

$MSAS_WORK_MDL/uucurr_pre_000000_002500_lo2  52  5x18  41_1 99  9080  900_00000000_modlerr 

#  rep 

$MSAS_WORK_MDL/ vveurr_pre_0 00 00 0_002 5 00_lo2  52  5x1 8  4 1_1 9 9  90 8 0  90 0_00 00 00 00_modlerr 

#  fi 

#  echo  " .  ready" 

FCST=${ifrcst} 

echo  'FORECAST  DATE  ’${FCST} 


print  -n  "  Checking  first  guess  fields  valid  tau  024" 

Cd  $ { RESTART_DIR} 

#  if  [  !  -r  seatmp_sf c_000000_000000_lo2  52  5x18  41_$ { FCST }_002  4  0000_f cstf Id  ]; 
then 

#  rep  $MSAS_WORK_MDL/hycom2dt_000000_$ { FCST }_0024 0000 .A 
seatmp_sf c_000000_000000_lo2525xl841_$ { FCST }_0024 0000_f cstf Id 

/ bi n / cp  $MSAS_LOCA_MDL/hycom2 dt_0 0  0  0  0  0_$ { FC  ST } _0  0  2  4  0  0  0  0 . A 
seatmp^sf c_000000_000000_lo2525xl841_$ { FCST }_0024 0000_f cstf Id 

#  fi 

#  if  [  !  -r  seahgt_sf c_000000_000000_lo2  52  5x18  41_$ { FCST }_002  4  0000_f cstf Id  ]; 
then 

#  rep  $MSAS_WORK_MDL/hycom2dSf sd_000000_$ { FCST }_0 02 4  0000 . A 
seahgt_sfc_000000_000000_lo2525xl841_$ {FCST }_00240000_fcstfld 

/bin/cp  $MSAS_LOCA_MDL/hycom2dSf sd_000000_$ {FCST }_0 02 40000 . A 
seahgt_sf c_000000_000000_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 
#for  tau  in  00030000  00060000  00090000  00120000  00150000  00180000  00210000 
00240000 
#do 

#  /bin/cp  $MSAS_LOCA_MDL/hycom2dSf sd_000000_$ { FCST }_$ { tau } .A 
seahgt_sf c_000000_000000_lo2525xl841_$ { FCST }_$ { tau}_f cstf Id 

#done 

#  fi 

#  if  [  !  -r  seaice_sf c_000000_000000_lo252 5x18 41_$ { FCST }_002 4  0000_f cstf Id  ]; 

then 

#  rep  $MSAS_WORK_MDL/hycom2dice_000000_$ { FCST }_0 02 4  0000 . A 
seaice_sfc_000000_000000_lo2525xl841_$ {FCST }_00240000_fcstfld 

/bin/cp  $MSAS_LOCA_MDL/hycom2 di ce_0 00000_${FCST}_00240000.A 
seaice_sf c_000000_000000_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 

#  fi 

#  if  [  !  -r  mixlyr_sf c_000000_000000_lo2 52 5x18 41_$ { FCST }_002 4  0000_f cstf Id  ]; 

then 

#  rep  $MSAS_WORK_MDL/hycom2dmix_000000_$ {FCST }_0 02 40000 . A 
mixlyr_sf c_000000_000000_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 

/bin/cp  $MSAS_LOCA_MDL/hycom2dmix_000000_$ {FCST }_0024 0000 .A 
mixlyr_sf c_000000_000000_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 

#  fi 

#  if  [  !  -r  seahgt_sf c_000000_000000_lo2 52 5x18 41_$ { FCST }_00000000_meanf Id  ]; 

then 

/bin/cp  $MSAS_WORK_MDL/RS_0 94_$ { T }_$ { REG1 } 0 . 08 . A 
seahgt_sf c_000000_000000_lo2525xl841_$ { FCST }_00 0000 00_meanf Id 

#  rep  $MSAS_WORK_MDL/ 056_archMNA. 0009_0013_fsd. A 

seahgt_sf c_000000_000000_lo2525xl841_$ { FCST }_00 0000 00_meanf Id 

#  rep  $MSAS_WORK_MDL/micomecmwf . year04_05_$ { REG } . A 

seahgt_sf c_000000_000000_lo2525xl841_$ { FCST }_00 0000 00_meanf Id 

#  rep  $MSAS_WORK_MDL/ zero_0 08_hycom . A 

seahgt^sf c_000000_000000_lo2525xl841_$ { FCST }_00 0000 00_meanf Id 

#  fi 

#  if  [  !  -r  densty_pre_000000_002500_lo2 52 5x18 41_$ { FCST }_002 4  0000_f cstf Id  ]; 

then 

##  rep  $MSAS_WORK_MDL/hycom3dden_002500_$ {FCST }_0 02 40000 . A 

densty_pre_000000_002  500_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 

#  /bin/cp  $MSAS_LOCA_MDL/hycom3dden_002500_$ { FCST }_0 02 4  0000 . A 
densty_pre_000000_002  500_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 
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#  fi 

#  if  [  !  -r  seatmp_pre_000000_002500_lo2  52  5x18  41_$ { FCST }_002  4  0000_f cstf Id  ]; 
then 

#  rep  $MSAS_WORK_MDL/hycom3dt_002500_$ {FCST }_0024  0000 .A 

seatmp_pre_000000_002500_lo2525xl841_$ { FCST }_0024 0000_f cstf Id 

/bin/cp  $MSAS_LOCA_MDL/hycom3dt_002  500_$ {FCST }_0024 0000 .A 
seatmp_pre_000000_002500_lo2525xl841_$ { FCST }_00240000_fcstfld 

#  fi 

#  if  [  !  -r  salint_pre_000000_002500_lo2 52 5x18 41_$ { FCST }_002 4  0000_f cstf Id  ]; 

then 

#  rep  $MSAS_WORK_MDL/hycom3ds_002500_$ { FCST }_00240000 .A 

salint_pre_000000_002  500_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 

/bin/cp  $MSAS_LOCA_MDL/hycom3ds_002 500_$ {FCST }_0024  0000 .A 
salint_pre_000000_002  500_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 

#  fi 

#  if  [  !  -r  uucurr_pre_000000_002500_lo2  52  5x18  41_$ { FCST }_002 4  0000_f cstf Id  ]; 

then 

#  rep  $MSAS_WORK_MDL/hycom3du_002  500_$ {FCST }_0024 0000 .A 

uucurr_pre_000000_002500_lo2525xl841_$ { FCST }_0024 0000_f cstf Id 

/bin/cp  $MSAS_LOCA_MDL/hycom3du_002 500_$ {FCST }_0024 0000 .A 
uucurr_pre_000000_002500_lo2525xl841_$ { FCST }_0024 0000_f cstf Id 

#  fi 

#  if  [  !  -r  vvcurr_pre_000000_002500_lo2 52 5x18 41_$ { FCST }_002 4  0000_f cstf Id  ]; 

then 

#  rep  $MSAS_WORK_MDL/hycom3dv_002500_${FCST}_00240000.A 

vvcurr_pre_000000_002  500_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 

/bin/cp  $MSAS_LOCA_MDL/hycom3dv_002 500_$ {FCST }_0024 0000 .A 
vvcurr_pre_000000_002  500_lo2525xl 84 1_$ { FCST }_0024  0000_f cstf Id 

#  fi 

if  [  !  -r  seahgt_sf c_000000_000000_lo2525xl 84 1_$ { FCST }_0024 0000_f cstf Id  ]; 

then 

echo  "seahgt  sfc  forecast  does  not  exist  for  "${FCST} 

exit 

fi 

if  [  !  -r  seatmp_sf c_000000_000000_lo2525xl 84 1_$ { FCST }_0024 0000_f cstf Id  ]; 

then 

echo  "seatmp  sfc  forecast  does  not  exist  for  "${FCST} 

exit 

fi 

if  [  !  -r  seatmp_pre_000000_002 500_lo2525xl 84 1_$ { FCST }_0024 0000_f cstf Id  ]; 

then 

echo  "seatmp  forecast  does  not  exist  for  "${FCST} 

exit 

fi 

if  [  !  -r  salint_pre_000000_002 500_lo2525xl 84 1_$ {FCST }_0024 0000_f cstf Id  ]; 

then 

echo  "salint  forecast  does  not  exist  for  "${FCST} 

exit 

fi 

if  [  !  -r  uucurr_pre_000000_002 500_lo2525xl 84 1_$ { FCST }_0024 0000_f cstf Id  ]; 

then 

echo  "uucurr  forecast  does  not  exist  for  "${FCST} 

exit 

fi 

if  [  !  -r  vvcurr_pre_000000_002 500_lo2525xl 84 1_$ { FCST }_0024 0000_f cstf Id  ]; 

then 

echo  "vveurr  forecast  does  not  exist  for  "${FCST} 


exit 

fi 

echo  " .  ready" 

print  -n  "  Checking  binary  files  " 
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if  [  !  -r  $BIN_DIR/ncoda_prep  ] ;  then 
echo  "Error:  ncoda_prep  not  found" 
fi 

if  [  !  -r  $BIN_DIR/ncoda  ] ;  then 
echo  "Error:  ncoda  not  found" 
fi 

if  [  !  -r  $BIN_DIR/ncoda_post  ] ;  then 
echo  "Error:  ncoda_post  not  found" 
fi 

echo  " .  ready" 

echo  "  " 
echo  "  " 

echo  "================================: 

echo  "  " 

echo  "  Data  preparation  complete:" 
echo  "  " 

echo  "================================: 


function  setup_processors  { 

# 

#  Set  OpenMP  environment 

#  Set  number  of  OpenMP  processors 

# 

#  Note:  The  pre  and  post  analysis  use  OpenMP  and  the 

#  analysis  uses  MPI  only. 

# 

#  Set  number  of  MPI  processors 

#  MPICMD= 'mpirun . lsf  ' 

#  MPICMD='poe  ' 
nprc=4 1 6 

#  Set  number  of  OpenMP  processors 

mprc=l 

# 

#  Set  OpenMP  environment 

# 

( (SLAVE=512*1024*1024 ) ) 
export  OMP_NUM_THREADS=$ {mprc } 

#  IBM 

export  XLSMPOPTS="stack=$SLAVE" 

#  XT5 

#  export  MPSTKZ=512M 

export  OMP_SCHEDULE="DYNAMIC, 1 " 

# 

#  export  OMPCMD= ' timex  -p  -mt  ' 

#  export  MPICMD= ' timex  -p  -mt  'poe 

#  export  MPICMD= ' timex  -p  -mt  'mpirun. lsf 

#  XT5 

export  OMPCMD= ' aprun  -n  1  -d  ' $ { OMP_NUM_THREADS } 
export  MPICMD= ' aprun  -n  '${nprc} 


echo  "  " 

echo  "##  PE  Environment  Overview  {" 
printenv  |  grep  -e  ,AOMP_'  -e  ' AMP '  -e 
'  AMPI  ' 


echo  "##  }" 
echo  "" 

# 

#  Print  the  hard  and  soft  limits. 

# 


DSM '  -e  ' A CHUNK '  - 


'  A  PAGE '  -e 
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ulimit  -aH 
ulimit  -aS 


function  setup_namelists  { 

# 

#  Setup  the  namelist  files 

# 

rm  -f  odsetnl 
rm  -f  gridnl 
rm  -f  oanl 

# 

#  data  path  settings 

# 

.  /u/home/$ {USER} /$ {REGO } /ncoda/$ { EXPT } /$ { REG1 }_$ {E } . odsetnl 

# 

#  grid  definition  namelist  settings 

#  kko  is  the  number  of  vertical  levels;  levels  are  defined  in  oanl 

# 

.  / u/home/ $ {USER} /$ { REGO } /ncoda/$ { EXPT } /$ { REG1 }_$ { E } .gridnl 

# 

#  ocean  analysis  namelist  settings 

#  for  CH  assimilation  of  altimeter  SSH,  set  direct  =  .true. 

#  for  MODAS  assimilation  of  altimeter  SSH,  set  modas  =  .true. 

#  for  SST  assimilation,  set  st_asm  =  .true. 

#  for  potential  temperature  analysis,  set  pt_anl  =  .true. 

#  (assumes  potential  temperature  background  fields) 

# 

.  / u/home/ $ {USER} /$ { REGO } /ncoda/$ { EXPT } /$ { REG1 }_$ { E } .oanl 

# 

echo  "  " 

} 


function  run_analysis  { 
typeset  analysis_dtg 
typeset  -i  cycle^index 

analysis_dtg=$ { idtg} 
cycle_index=$ { inumd} 

# 

#  Remove  the  work  files  created  by  the  analysis 

# 

/bin/rm  -f  pout* 


echo  "  Start  Ocean  Prep  2D  at  "  'date  +%T'  »  poutl 
echo  "  "  »  poutl 

${OMPCMD}  $BIN_DIR/ncoda  prep  2D  gridnl  $analysis  dtg  >>  poutl 


if  [  $?  -ne 
echo 
echo 
echo 
echo 
echo 
echo 
fatal 
else 
echo 
echo 
echo 
echo 
echo 
echo 
fi 


0  ] ;  then 
»  poutl 
»  poutl 

Error;  Ocean  Prep  2D  failed  at 
Check  poutl"  »  poutl 
»  poutl 
»  poutl 


»  poutl 
»  poutl 
End  Ocean 
Completed 
»  poutl 
»  poutl 


Prep  2D  at  " 
successfully" 


'date  +%T '  »  poutl 


date  +%T ' 
»  poutl 


>>  poutl 
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echo  "  Start  Ocean  Analysis  2D  at  "  'date  +%T'  »  pout2 
echo  "  "  >>  pout2 

${MPICMD}  $BIN_DIR/ncoda  2D  gridnl  $analysis_dtg  »  pout2 
if  [  $?  -ne  0  ] ;  then 
echo  "  "  »  pout2 
echo  "  "  »  pout2 

echo  "  Error:  Ocean  Analysis  2D  failed  at  "  'date  +%T'  »  pout2 
echo  "  Check  pout2"  »  pout2 
echo  "  "  »  pout2 
echo  "  "  »  pout2 
fatal 
else 

echo  "  "  »  pout2 
echo  "  "  »  pout2 

echo  "  End  Ocean  Analysis  2D  at  "  'date  +%T'  >>  pout2 

echo  "  Completed  successfully"  »  pout2 

echo  "  "  »  pout2 
echo  "  "  »  pout2 


echo  "  Start  Ocean  Post  2D  at  "  'date  +%T'  »  pout3 
echo  "  "  >>  pout3 

${OMPCMD}  $BIN_DIR/ncoda_post  2D  gridnl  $analysis_dtg  >>  pout3 
if  [  $?  -ne  0  ] ;  then 
echo  "  "  »  pout3 
echo  "  "  »  pout3 

echo  "  Error:  Ocean  Post  2D  failed  at  "  'date  +%T'  >>  pout3 
echo  "  Check  pout3"  »  pout3 
echo  "  "  »  pout3 
echo  "  "  »  pout3 
fatal 
else 

echo  "  "  »  pout3 
echo  "  "  »  pout3 

echo  "  End  Ocean  Post  2D  at  "  'date  +%T'  >>  pout3 
echo  "  Completed  successfully"  >>  pout3 
echo  "  "  »  pout3 
echo  "  "  »  pout3 


echo  "  Start  Ocean  Prep  3D  at  "  'date  +%T'  »  pout4 
echo  "  "  »  pout4 

${OMPCMD}  $BIN_DIR/ncoda_prep  3D  gridnl  $analysis_dtg  >>  pout4 
if  [  $?  -ne  0  ] ;  then 
echo  "  "  »  pout4 
echo  "  "  »  pout4 

echo  "  Error:  Ocean  Prep  3D  failed  at  "  'date  +%T'  >>  pout4 

echo  "  Check  pout4"  »  pout4 

echo  "  "  »  pout4 
echo  "  "  »  pout4 
fatal 
else 

echo  "  "  »  pout4 
echo  "  "  »  pout4 

echo  "  End  Ocean  Prep  3D  at  "  'date  +%T'  >>  pout4 
echo  "  Completed  successfully"  >>  pout4 
echo  "  "  »  pout4 
echo  "  "  »  pout4 


echo  "  Start  Ocean  Analysis  3D  at  "  'date  +%T'  »  pout5 
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echo  "  "  >>  pout5 

${MPICMD}  $BIN_DIR/ncoda  3D  gridnl  $analysis_dtg  »  pout5 
if  [  $?  -ne  0  ] ;  then 

if  ( (  $cycle_index  ==  1  ) ) ;  then 
echo  "  "  »  pout5 
echo  "  "  »  pout5 

echo  "  Error:  Ocean  Analysis  3D  failed  at  "  'date  +%T'  >>  pout5 

echo  "  Check  pout5"  »  pout5 

echo  "  This  is  not  a  fatal  condition,  continuing  ..."  »  pout5 
echo  "  "  »  pout5 
echo  "  " 
else 

echo  "  "  »  pout5 
echo  "  "  »  pout5 

echo  "  Error:  Ocean  Analysis  3D  failed  at  "  'date  +%T'  >>  pout5 

echo  "  Check  pout5"  »  pout5 

echo  "  "  »  pout5 
echo  "  "  »  pout5 
fatal 
fi 

else 

echo  "  "  »  pout5 
echo  "  "  »  pout5 

echo  "  End  Ocean  Analysis  3D  at  "  'date  +%T'  >>  pout5 

echo  "  Completed  successfully"  »  pout5 

echo  "  "  »  pout5 
echo  "  "  »  pout5 


echo  "  Start  Ocean  Post  3D  at  "  'date  +%T'  »  pout6 
echo  "  "  >>  pout 6 

${OMPCMD}  $BIN_DIR/ncoda_post  3D  gridnl  $analysis_dtg  >>  pout6 
if  [  $?  -ne  0  ] ;  then 

if  ( (  $cycle_index  ==  1  ) ) ;  then 
echo  "  "  »  pout6 
echo  "  "  »  pout6 

echo  "  Error:  Ocean  Post  3D  failed  at  "  'date  +%T'  >>  pout6 

echo  "  Check  pout6"  »  pout6 

echo  "  This  is  not  a  fatal  condition,  continuing  ..."  »  pout6 
echo  "  "  »  pout6 
echo  "  "  »  pout6 
else 

echo  "  "  »  pout6 
echo  "  "  »  pout6 

echo  "  Error:  Ocean  Post  3D  failed  at  "  'date  +%T'  >>  pout6 

echo  "  Check  pout6"  »  pout6 

echo  "  "  »  pout6 
echo  "  "  »  pout6 
fatal 
fi 

else 

echo  "  "  »  pout6 
echo  "  "  »  pout6 

echo  "  End  Ocean  Post  3D  at  "  'date  +%T'  >>  pout6 
echo  "  Completed  successfully"  »  pout6 
echo  ""  >>  pout 6 
echo  ""  >>  pout6 
fi 

# 

#  rename  diagnostic  files  (creation  controlled  by  debug  options  in  oanl) 

# 

mv  fort . 25  $ { SAVE_DIR} /$ { REG1 } _ $ { E } _ $ { idtgtod}_$ {analysis_dtg} . rpr 

mv  fort . 27  $ { SAVE_DIR} / $ { REG1 } _ $ { E } _ $ { idtgtod}_$ {analysis_dtg} . vfy 
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# 

# 

# 


mv 

fort . 

31 

mv 

fort . 

32 

mv 

fort . 

33 

mv 

fort . 

34 

mv 

fort . 

35 

mv 

fort . 

36 

mv 

fort . 

37 

mv 

fort . 

38 

mv 

fort . 

39 

$ { SAVE_DIR} /$ { REG1 } _ $ { E } _ $ { idtgtod }_$ {analysis_dtg} . syn 

$ { SAVE_DIR} / $ { REG1 } _ $ { E }_$ { idtgtod} _$ { analysi s_dtg } . rej 

$ { SAVE_DIR} /$ {REG1 } _ $ { E } _ $ { idtgtod }_$ { analysi s_dtg} .prf 

$ {SAVE_DIR} /$ { REG1 } _$ { E }_$ { idtgtod} _$ { analysi s_dtg } .gpt 

${SAVE_DIR} / $ { REG1 } _ $ { E }_$ { idtgtod }_$ { analysi s_dtg} .err 

$ { SAVE_DIR} /$ { REG1 } _ $ { E } _ $ { idtgtod }_$ { analysi s_dtg} .mvo 

$ { SAVE_DIR} /$ {REG1 } _ $ { E } _ $ { idtgtod }_$ { analysi s_dtg} . drc 

${SAVE_DIR} /${REG1 }_$ { E }  $ { idtgtod} _$ {analysi s_dtg} . lyp 
${SAVE_DIR} / $ { REG1 }_$ { E }_$ { idtgtod }_$ { analysi s_dtg} . fix 


combine  program  outputs  into  single  file  and  finish 


/bin/rm  $analysis_dtg . out 

cat  poutl  pout2  pout3  pout4  pout5  pout6  >>  $analysis_dtg . out 

echo  " -  odsetnl  - "  »  $analysis_dtg . out 

cat  odsetnl  »  $analysis_dtg. out 

echo  " -  gridnl  - "  >>  $analysis_dtg. out 

cat  gridnl  »  $analysis_dtg . out 

echo  " - oanl - "  >>  $analysis_dtg.  out 

cat  oanl  »  $analysis_dtg . out 

echo  "  Analysis  ended  at  "  'date  +%T' 

echo  "  Output  file:  $analysis_dtg . out  in  "  'pwd' 

echo  "  "  'date  -u  +'%Y%m%d'00' 

echo  "  " 

echo  "  " 

echo  "====================================================" 

mv  $analysis_dtg. out  $ { SAVE_DIR} / $ { REG1 } $ { E } $ { idtgtod}_$ { analysis_dtg } . out 


function  cleanup  { 

cd  $RUN_DIR/analysis 

#  /bin/rm  -f  coda.*obs& 

#  /bin/ls  |  egrep  ' AMVOI '  |  xargs  /bin/rm 

/bin/rm  MVOI*& 

#  rm  -f  pout* 

} 


# 

# - 

# - 

# 

#  Begin  script  execution,  check  command  line  arguments 

# 

# - 

# - 

# 

#if  [  "$#"  -ne  3  ];  then 

#  clear 

#  usage 

#  print  ERROR:  incorrect  number  of  input  args 

#  exit  1 
#else 

FIRST_DTG="$ { idtg} " 

NUM_DAYS=" $ { inumd} " 

FCST="${ifrcst}" 

#fi 

if  [[  $ { FIRST_DTG }  !=  +([0-9])  ]];  then 
clear 
usage 

print  ERROR:  arg  1  is  not  an  integer 
exit  1 
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fi 

if  [[  ${NUM_DAYS]  !=  +([0-9])  ]];  then 
clear 
usage 

print  ERROR:  arg  2  is  not  an  integer 
exit  1 
fi 

if  [[  $ { FCST }  ! =  +([0-9])  ]];  then 
clear 
usage 

print  ERROR:  arg  3  is  not  an  integer 
exit  1 
fi 

if  [  $ { #FIRST_DTG }  -ne  10  ];  then 
clear 
usage 

print  ERROR:  incorrect  format  for  arg  1 
exit  1 
fi 

if  [  $ { #FCST }  -ne  10  ];  then 
clear 
usage 

print  ERROR:  incorrect  format  for  arg  3 
exit  1 
fi 

# 

#  The  batch  process  may  set  file  creation  mask  to  600,  change 

#  this  with  the  mask. 

# 

umask  022 

# 

#  Define  some  useful  variables  for  the  model  name,  dimensions, 

#  and  number  of  hours  in  update  cycle. 

# 

MODEL=hycom 
UPDATE_CYCLE=24 
setup_paths  $ [MODEL] 

# 

#  Setup  date-time-groups  for  the  script 

# 

let  "day_index=$ { NUM_DAYS }  -  1" 
export  CRDATE=$ [FIRST_DTG: -19 99060600} 

prev_dtg= ' $DDTG  -d  -1  2>  /dev/null' 

#  prev_dtg= ' $DDTG  -d  -7  2>  /dev/null' 

LAS T_DTG= ' $ DDTG  -d  ${day_index}  2>  /dev/null' 
curr_dtg= ' $DDTG  2>  /dev/null' 


echo 

=  =  =  = 

— 

— 

echo 

II 

II 

echo 

II 

II 

echo 

"CODA 

analysis  processor (ncoda 

$ [REG1 } .job) " 

echo 

II 

The  current  date/time  is" 

$ (  date) 

echo 

II 

II 

echo 

II 

II 

echo 

II 

Dates :  " 

echo 

II 

First  cycle  date 

:  $FIRST  DTG" 

echo 

II 

Last  cycle  date 

:  $LAST  DTG" 
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echo  "  Forecast  date  :  $FCST" 

echo  "  " 
echo  "  " 

echo  "  Path  setting:" 

echo  "  Binary  locations:  $BIN_DIR" 

echo  "  Scratch  directory:  $SCRATCH_DIR" 

#  echo  "  Model  directory:  $MODEL_DIR" 

echo  "  Restart  directory:  $RESTART_DIR" 

echo  "  CODA  database  directory:  $DB_DIR  " 

echo  "  CODA  obs  data  directory:  $OBS_DIR  " 

echo  "  " 
echo  "  " 

echo  "====================================================" 


setup_dirs 

# 

#  Clean  the  temporary  directory  structure 

# 

if  [[  !  -e  $RUN_DIR  ]];  then 

echo  "The  run  directory  is  missing,  can  not  proceed:  $RUN_DIR" 
exit  1 
fi 

if  [ [  !  -e  $RUN_DIR/restart  ] ] ;  then 

echo  "The  Restart  directory  is  missing,  can  not  proceed:  $RUN_DIR/ restart" 
exit  1 
fi 

cd  $RUN_DIR 

setup_data  $ { FIRST_DTG}  ${NUM_DAYS}  ${FCST} 
setup_processors 

# 

#  Execute  the  coda  mvoi  programs 

# 

echo  "  " 
echo  "  " 

echo  "Start  Analysis" 

echo  "  Directories  created" 

echo  "  Date-time-group :  $curr_dtg" 

echo  "  MPI  processes:  $nprc" 

echo  "  OpenMP  threads:  $mprc" 

echo  "  Initial  directory  "  'pwd' 

echo  "  " 

cd  $RUN_DIR/analysis 
setup_namelists 

# 

#  loop  over  update  cycles 

# 

i=l 

h=0 

num_cycles=$ { NUM_DAYS } 

#while  ( (  $i  <=  $num_cycles  ) )  do 

# 

#  set  analysis  dtg 

# 

curr_dtg= ' $DDTG  -h  $h  2>  /dev/null' 
echo  'UPDATE  CYCLE  :  ' $UPDATE  CYCLE 
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echo  'curr_dtg  :  '${curr_dtg}  $ { i } 
idtg=${idtg} 

#  execute  ncoda  mpi  mvoi  programs 

run_analysis  ${idtg}  ${i} 

#  run_analysis  ${curr_dtg}  $ { i } 

#  update  counters 

#  (  (  i  =  i  +  1  )  ) 

#  (  (  h  =  h  +  $UPDATE_CYCLE  )  ) 

#done 

♦cleanup 

# 

cd  /scr/$ {USER} /hycom/ $ { REGO } /$ { EXPT } /logs 
touch  done_$ { REG1 }_$ { idtgtod}_$ { curr_dtg} 

/u/home/wallcraf /bin/q_navo  $ { E } ncoda2arch_$ { idtgtod}_$ { curr_dtg} . csh 

# 

/u/home/wallcraf /bin/ q_navo  $ {E }plot_$ { REG1 }_$ { idtgtod}_$ { curr_dtg} . csh 
echo  "" 

echo  "====================================================" 

exit  0 

#!/bin/csh  -f 
#PBS  -N  999pbs 
#PBS  -j  oe 

#PBS  -o  999pbsXX.log 
#PBS  -1  mppwidth=l 
#PBS  -1  mppnppn=l 
#PBS  -1  mppmem=16gb 
#PBS  -1  walltime=2 : 00 : 00 
#PBS  -W  umask=027 
#PBS  -A  NRLSS018 
#PBS  -q  internal3d 
# 

#BSUB  -J  727ncoda2arch_all 

#BSUB  -n  1 

#BSUB  -W  2:30 

#BSUB  -R  "span [ptile=l] " 

#BSUB  -q  internal3d 
#BSUB  -P  NAVOSHYC 
# 

set  echo 
set  time  =  1 

set  BIN=/u/home/$ {user } /bin 
set  pget=~wallcraf /bin/pget 
# 

set  OS='uname' 
switch  ($OS) 
case  1  Linux  1 : 

which  aprun 

if  (!  $status)  then 

set  APRUN=' aprun  -n  1  -m  16g  ' 
set  SRC=/ u/home/wallcraf /hycom/ALLcnl 
else 

set  APRUN= ' ' 

set  SRC=/ u/home/wallcraf /hycom/ ALL 
endif 
breaksw 
case  ' AIX ' 
set  APRUN= ' ' 

set  SRC= /u/home/wallcraf /hycom/ ALL 
breaksw 
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default : 
set  APRUN= ' ' 

set  SRC= /u /home /waller a f /hycom/ ALL 
endsw 
# 

C 

C  -  Remap  an  archive  file  to  an  NCODA  analysis,  new  layer  depths. 

C 

setenv  pget  cp 
setenv  pput  cp 

setenv  pget  -wallcraf /bin/pget 
setenv  pput  -wallcraf /bin/pput 
C 


C - E 

is 

experiment  number 

C - REG 

is 

region  identifier 

C - T 

is 

topog.  identifier 

C - EN 

is 

ncoda  experiment  number 

C - REGN 

C 

is 

ncoda  region  identifier 

# 

#  analysis 

hour 

# 

setenv  HR  18 

# 

setenv  idtg  20080501 
setenv  idtgl  ${idtg} 
setenv  idtgtod  20080501 
setenv  idtgtodpl  20080501 
setenv  idtgtodml  20080501 

setenv  idtgtodplt  ' ${BIN} /addndays  yyyymmdd  ${idtgtod}  -1' 
#setenv  idtgmax  '$ {BIN} /addndays  yyyymmdd  ${idtg}  +0' 

# 

setenv  E  727 
set  REG=GLBa0 . 08 
set  EXPT=expt_72 . 7 
set  case=nowcast 
setenv  T  09 

setenv  MVOI  /u/home/$ {user } /hycom/$ {REG} /${ EXPT } /mvoi 
set  TOPO=/ u/home/ $ {user } /hycom/ $ {REG} /topo 
set  RUN=/ u/home/ $ { user } /hycom/ $ { REG} /$ { EXPT } /mvoi 
set  LOGS=/scr/$ { user } /hycom/ $ { REG} / $ {EXPT } /logs 
set  INP=/ scr/$ {user } /hycom/ $ {REG} / $ { EXPT } /$ { case } 

# 

set  arcl=archv 
set  arc2=archv_l 
# 

cd  $ { INP } 

C 

touch  ncoda_archv_vel 
/bin/rm  ncoda_archv_vel 
touch  ncoda_archv 
/bin/rm  ncoda_archv 
# 

foreach  REGN (  ANTarc  MERatl  MERpal  MERinl  ARCatl  ARCpac  ARCocn) 

# 

START : 

if (  -e  $ { LOGS } / done_$ {REGN }_$ { idtgtodpl } $ { HR}_$ { idtg} $ { HR} )  then 
goto  RUN 
else 

echo  'Sleeping' 
sleep  60 
goto  START 
endif 
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RUN: 

# 

if  (${ 

REGN }  ==  "ARCocn 

set 

SI ZN=" 1600x1297" 

set 

ilstn=2459 

set 

j lstn=2650 

set 

idmn=l 600 

set 

jdmn=1297 

endif 

if  (${ 

REGN}  ==  "ARCatl 

set 

SI ZN=" 1490x0551" 

set 

ilstn=2440 

set 

j lstn=2150 

set 

idmn=1490 

set 

jdmn=551 

endif 

if  (${ 

REGN}  ==  "ARCpac 

set 

SI ZN=" 1335x0 551" 

set 

ilstn=730 

set 

j lstn=2150 

set 

idmn=1335 

set 

jdmn=551 

endif 

if  (${ 

REGN}  ==  "ANTarc 

set 

SI ZN=" 3827x0411" 

set 

il stn=82  4 

set 

j lstn=l 

set 

idmn=3827 

set 

jdmn=411 

endif 

if  (${ 

REGN}  ==  "MERla" 

set 

SI ZN=" 4501x0511" 

set 

ilstn=3574 

set 

j lstn=360 

set 

idmn=4501 

set 

jdmn=511 

endif 

if  (${ 

REGN}  ==  "MER2b" 

set 

SI ZN=" 4501x0511" 

set 

ilstn=3574 

set 

j lstn=820 

set 

idmn=4501 

set 

jdmn=511 

endif 

if  (${ 

REGN}  ==  "MER3c" 

set 

SI ZN=" 4501x0511" 

set 

ilstn=3574 

set 

j lstn=1280 

set 

idmn=4501 

set 

jdmn=511 

endif 

if  (${ 

REGN}  ==  "MER4d" 

set 

SI ZN=" 4501x0461" 

set 

ilstn=3574 

set 

j 1 stn=17  40 

set 

idmn=4501 

set 

j  dmn=4  61 

endif 

if  (${ 

REGN}  ==  "MERatl 

set 

SI ZN=" 1751x1841" 

set 

ilstn=2349 

set 

j lstn=360 

set 

idmn=1751 

)  then 


)  then 


)  then 


)  then 


then 


then 


then 


then 


)  then 
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set  jdmn=1841 
endif 
# 

if ( $ { REGN }  ==  "MERpal")  then 
set  SI ZN=" 2525x1841" 
set  ilstn=199 
set  jlstn=360 
set  idmn=2525 
set  jdmn=1841 
endif 
# 

if ($ {REGN}  ==  "MERinl")  then 
set  SI ZN=" 1313x1569" 
set  ilstn=3949 
set  jlstn=360 
set  idmn=1313 
set  jdmn=1569 
endif 
# 

setenv  EN  727 

set  EXPTN=expt_72 . 7 

setenv  TN  ${T} 

set  T0P0N=/u/home/$ {user} /hycom/ $ {REG} / subtopo/ $ { REGN } 

# 

set  NCODA=/scr/$ {user} /${REGN}0 . 08/${EXPTN} /restart 
mkdir  -p  ${ INP} /.. /data/incup 
# 

# 

set  typet=seatmp_pre_000000_002500_lo$ { SIZN } 
set  types=salint_pre_000000_002500_lo$ { SIZN } 
set  typeu=uucurr_pre_000000_002500_lo$ { SIZN } 
set  typev=vvcurr_pre_000000_002500_lo$ { SIZN } 
set  typep=lyrprs_pre_000000_002500_lo$ { SIZN } 
set  typei=seaice_sf c_000000_000000_lo$ { SIZN } 

C 

touch  zi42.txt 

if  ( —  z  zi42.txt)  then 

#  ${pget}  $ { RUN } /zi42 . txt  zi42.txt  & 

/bin/cp  $ { RUN } /zi42 . txt  zi42.txt  & 

endif 

C 

touch  regional . depth . a  regional . depth . b 
if  ( —  z  regional . depth . b)  then 

#  ${pget}  $ { TOPO} /depth_$ { REG}_$ {T } . b  regional . depth . b  & 
cp  $ { TOPO} /depth_$ { REG}_$ {T } . b  regional . depth . b  & 

endif 

if  ( —  z  regional . depth . a)  then 

#  ${pget}  ${TOPO}/depth_${REG}_${T} .a  regional . depth . a  & 
cp  $ { TOPO} /depth_$ { REG}_$ {T } . a  regional . depth . a  & 

endif 

C 

touch  regional . grid . a  regional . grid . b 
if  ( —  z  regional . grid. b)  then 

#  ${pget}  ${ TOPO} /regional . grid. b  regional . grid. b  & 
cp  ${ TOPO } /regional . grid. b  regional . grid. b  & 

endif 

if  ( —  z  regional . grid. a)  then 

#  ${pget}  ${ TOPO} /regional . grid. a  regional . grid. a  & 
cp  ${ TOPO } /regional . grid. a  regional . grid. a  & 

endif 

touch  iso. sigma. a  iso. sigma. b 
if  ( —  z  iso. sigma. b)  then 

cp  ${ INP} /.. /data/iso . sigma .b  iso. sigma. b  & 
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endif 

if  ( - z  iso. sigma. a)  then 

cp  ${ INP} /.. /data/iso . sigma . a  iso. sigma. a  & 
endif 
# 

#  subregion  depth  files 

# 

touch  subregional . depth . $ { REGN } . a 

if  ( - z  subregional . depth .${ REGN }. a)  then 

#  ${pget}  $ { TOPON } / depth_$ { REGN } 0 . 0 8_$ { TN } . a  subregional . depth .${ REGN }. a  & 
/bin/cp  $ {TOPON} /depth_$ { REGN} 0 . 08_$ { TN} . a  subregional . depth .${ REGN }. a  & 

endif 

C 

touch  ncoda_archv_vel 

touch  ncoda_archv 

#if  ( —  z  ncoda_archv_vel )  then 

#  cp  -wallcraf /hycom/ALL/archive/src/ncoda_archv_vel  .  & 

cp  $ { SRC } /archive/src/ncoda_archv  .  & 

#endif 

wait 

chmod  a+rx  ncoda_archv_vel 
chmod  a+rx  ncoda_archv 
C 

# 

/bin/rm  $ { arcl } . $ { idtgl } . a  $ { arcl } . $ { idtgl } . b 

touch  $ { arcl } . $ { idtgtodml } $ { HR}_$ { idtgl } $ { HR} .a 

if  ( - z  ${ arcl }.${ idtgtodml }${ HR}_$ { idtgl }${ HR} . a)  then 

echo  ${ arcl }.${ idtgtodml }${ HR}_$ { idtgl }${ HR} . a  is  not  there 
else 

In  -s  $ {arcl }.${ idtgtodml }$ {HR}_$ { idtgl }$ {HR} . a  ${ arcl }.${ idtgl }. a 
endif 

touch  $ { arcl } . $ { idtgtodml } $ { HR}_$ { idtgl } $ { HR} . b 

if  ( - z  ${ arcl }.${ idtgtodml }${ HR}_$ { idtgl }${ HR} . b)  then 

echo  ${ arcl }.${ idtgtodml }${ HR}_$ { idtgl }${ HR} . b  is  not  there 
else 

In  -s  $ {arcl }.${ idtgtodml }$ {HR}_$ { idtgl }$ {HR} .b  ${ arcl }.${ idtgl }. b 
endif 
C 

/bin/ rm  seatmp_$ { idtg }_ncoda_$ { REGN } 
touch  seatmp_$ { idtg}_ncoda_$ { REGN } 
if  ( - z  seatmp_$ { idtg}_ncoda_$ {REGN} )  then 
/bin/rm  seatmp_$ { idtg }_ncoda_$ { REGN } 

In  -s  $ {NCODA} /$ { typet }_$ { idtg} $ { HR}_00000000_analf Id 
seatmp_${ idtg  }_jncoda_${  REGN}  & 
endif 

/bin/rm  salint_$ { idtg }_ncoda_$ { REGN } 
touch  salint_$ { idtg }_ncoda_$ { REGN } 
if  ( - z  salint_$ { idtg }_ncoda_$ { REGN } )  then 
/bin/rm  salint_$ { idtg }_ncoda_$ { REGN } 

In  -s  $ {NCODA} /${ types }_${ idtg} ${ HR}_00000000_analf Id 
salint_$ { idtg }_ncoda_$ { REGN}  & 
endif 

/bin/rm  uucuri_$ { idtg }_ncoda_$ { REGN } 
touch  uucuri_$ { idtg }_ncoda_$ { REGN } 
if  ( - z  uucuri_$ { idtg }_ncoda_${ REGN} )  then 
/bin/rm  uucur i_$ { idtg }_ncoda_$ { REGN } 

In  -s  $ {NCODA} /${ typeu}_$ { idtg} ${ HR}_00000000_analinc 
uucuri_${  idtg }__ncoda_${ REGN}  & 
endif 

/bin/ rm  vvcur i_$ { idtg }_ncoda_$ { REGN } 
touch  vvcuri_$ { idtg }_ncoda_$ { REGN } 
if  ( —  z  vvcuri_$ { idtg }_ncoda_$ { REGN } )  then 
/bin/rm  vvcur i_$ { idtg }_ncoda_$ { REGN } 
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In  -s  $ {NCODA} /$ { typev}_$ { idtg} $ { HR}_00000000_analinc 
vvcuri_${ idtg }_ncoda_${ REGN}  & 
endif 

/bin/ rm  layprs_$ { idtg }_ncoda_$ { REGN } 
touch  layprs_$ { idtg }_ncoda_$ { REGN } 

if  (-z  layprs_$ { idtg }_ncoda_$ {REGN } )  then 
/bin/rm  laypr s_$ { idtg }_ncoda_$ { REGN } 

In  -s  $ {NCODA} /${ typep}_$ { idtg} ${ HR}_00000000_analinc 
layprs_$ { idtg }_ncoda_$ { REGN}  & 
endif 

/bin/ rm  seaice_$ { idtg }_ncoda_$ { REGN } 
touch  seaice_$ { idtg}_ncoda_$ { REGN } 
if  ( - z  seaice_$ { idtg }_ncoda_$ { REGN } )  then 
/bin/rm  seaice_$ { idtg }_ncoda_$ { REGN } 

In  -s  $ {NCODA} /${ typei }_${ idtg} ${ HR}_00000000_analf Id 
seaice_$ { idtg }_ncoda_$ { REGN}  & 
endif 


wait 


/bin/rm  -f  $ { arc2 } . $ { idtg } . a  $ { arc2 } . $ { idtg } . b 

# 

#  foreach  REGN 

end 

date 

C 

C  -  do  all  subregions  in  a  single  invokation  of  ncoda_archv_vel 

C 

date 

${APRUN}  . /ncoda_archv  «E-o-D 
$ {arcl } . $ { idtgl } . a 
${arc2} . $ { idtg} .a 

0  'intflg'  =  vertical  interpolation  flag  (0=T&S,  l=th&S) 

2  'isoflg'  =  preserve  isopycnal  layer  flag  ( 0=n , l=y, 2=y&layT, 3=y&isoT ) 

000  ' iexpt  '  =  experiment  number  xlO  (000=from  archive  file) 


3 

'yrflag'  = 

days  in  year  flag 

(0=360, 

1=366,  2=366 J1 , 

3= 

actual) 

4500 

'  idm 

=  longitudinal  array  size 

3298 

'  jdm 

=  latitudinal  array  size 

2460 

'itest  ' 

=  longitudinal  test  point 

(optional,  default 

0) 

1874 

'jtest  ' 

=  latitudinal  test  point 

(optional,  default 

0) 

32 

' kdm  '  = 

number  of  layers 

32 

' nhybrd ' 

=  number  of  hybrid 

levels 

(0=all  isopycnal) 

14 

' nsigma ' 

=  number  of  sigma 

levels 

(nhybrd-nsigma  z-levels) 

3.0 

' dpOO  ' 

=  deep  z-level 

spacing 

minimum  thickness 

(m) 

450.0 

' dpOOx  ' 

=  deep  z-level 

spacing 

maximum  thickness 

(m) 

1 . 18 

' dpOOf  ' 

=  deep  z-level 

spacing 

stretching  factor 

(1. 

0=const . 

0.5 

' dsOO  ' 

=  shallow  z-level 

spacing 

minimum  thickness 

(m) 

75.0 

'dsOOx  ' 

=  shallow  z-level 

spacing 

maximum  thickness 

(m) 

1 . 18 

' dsOOf  ' 

=  shallow  z-level 

spacing 

stretching  factor 

(1. 

0=const . 

1.0 

' dpOOi  ' 

=  deep  iso-pycnal 

spacing 

minimum  thickness 

(m) 

6.0 

' isotop ' 

=  shallowest  depth 

for  isopycnal  layers 

(m. 

<0  from 

0.03 

' deniso ' 

=  isopycnal  if  layer  is  within  deniso  of  target 

density 

34.0 

'thbase'  = 

reference  density 

(sigma 

units) 

1 

' vsigma ' 

=  spacially  varying 

isopycnal  target 

densities  (0 

28.10 

' sigma  ' 

=  layer 

1 

isopycnal 

target 

density 

(sigma 

units ) 

28.90 

' sigma  ' 

=  layer 

2 

isopycnal 

target 

density 

(sigma 

units) 

29.70 

' sigma  ' 

—  layer 

3 

isopycnal 

target 

density 

(sigma 

units ) 

30.50 

' sigma  ' 

=  layer 

4 

isopycnal 

target 

density 

(sigma 

units ) 

30.95 

' sigma  ' 

=  layer 

5 

isopycnal 

target 

density 

(sigma 

units) 

31 . 50 

' sigma  ' 

=  layer 

6 

isopycnal 

target 

density 

(sigma 

units ) 

32.05 

' sigma  ' 

=  layer 

7 

isopycnal 

target 

density 

(sigma 

units ) 

32.60 

' sigma  ' 

=  layer 

8 

isopycnal 

target 

density 

(sigma 

units ) 

33.15 

' sigma  ' 

=  layer 

9 

isopycnal 

target 

density 

(sigma 

units) 
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33.70 

' sigma  ' 

= 

layer 

10 

isopycnal 

target 

density 

(sigma 

units ) 

34.25 

' sigma  ' 

= 

layer 

11 

isopycnal 

target 

density 

(sigma 

units) 

34 . 75 

' sigma  ' 

= 

layer 

12 

isopycnal 

target 

density 

(sigma 

units ) 

35.15 

' sigma  ' 

= 

layer 

13 

isopycnal 

target 

density 

(sigma 

units) 

35.50 

' sigma  ' 

= 

layer 

14 

isopycnal 

target 

density 

(sigma 

units ) 

35.80 

' sigma  ' 

= 

layer 

15 

isopycnal 

target 

density 

(sigma 

units ) 

36.04 

' sigma  ' 

= 

layer 

16 

isopycnal 

target 

density 

(sigma 

units ) 

36.20 

' sigma  ' 

= 

layer 

17 

isopycnal 

target 

density 

(sigma 

units) 

36.38 

' sigma  ' 

= 

layer 

18 

isopycnal 

target 

density 

(sigma 

units ) 

36.52 

' sigma  ' 

= 

layer 

19 

isopycnal 

target 

density 

(sigma 

units ) 

36.62 

' sigma  ' 

= 

layer 

20 

isopycnal 

target 

density 

(sigma 

units ) 

36.70 

' sigma  ' 

= 

layer 

21 

isopycnal 

target 

density 

(sigma 

units) 

36.77 

' sigma  ' 

= 

layer 

22 

isopycnal 

target 

density 

(sigma 

units ) 

36.83 

' sigma  ' 

= 

layer 

23 

isopycnal 

target 

density 

(sigma 

units ) 

36.89 

' sigma  ' 

= 

layer 

24 

isopycnal 

target 

density 

(sigma 

units) 

36.97 

' sigma  ' 

= 

layer 

25 

isopycnal 

target 

density 

(sigma 

units ) 

37.02 

' sigma  ' 

= 

layer 

26 

isopycnal 

target 

density 

(sigma 

units) 

37.06 

' sigma  ' 

= 

layer 

27 

isopycnal 

target 

density 

(sigma 

units) 

37 . 10 

' sigma  ' 

= 

layer 

28 

isopycnal 

target 

density 

(sigma 

units ) 

37 . 17 

' sigma  ' 

= 

layer 

29 

isopycnal 

target 

density 

(sigma 

units) 

37.30 

' sigma  ' 

= 

layer 

30 

isopycnal 

target 

density 

(sigma 

units) 

37 . 42 

' sigma  ' 

= 

layer 

31 

isopycnal 

target 

density 

(sigma 

units ) 

37 . 48 

' sigma  ' 

= 

layer 

32 

isopycnal 

target 

density 

(sigma 

units ) 

0.5 

' hicemn ' 

= 

minimum 

ice  thickness  (m) 

seatmp_${ idtg}_ncoda_ANTarc 
salint_$ { idtg }_ncoda_ANTarc 
uucur i_$ { idtg }_ncoda_ANTarc 
vvcur i_$ { idtg }_ncoda_ANTarc 
layprs_${ idtg }_ncoda_ANTarc 
subregional . depth . ANTarc . a 
seaice_${ idtg}_ncoda_ANTarc 
zi42 . txt 

824  'ilstn  '  =  i-origin  of  ncoda  subregion 

1  'jlstn  '  =  j-origin  of  ncoda  subregion 

3827  ' idmn  '  =  i-extent  of  ncoda  subregion  (<=idm;  0  implies  idm) 

411  'jdmn  '  =  j-extent  of  ncoda  subregion  (<=jdm;  0  implies  jdm) 

42  ' kncoda '  =  number  of  ncoda  levels 

seatmp_$ { idtg }_ncoda_MERinl 
salint_$ { idtg }_ncoda_MERinl 
uucur i_$ { idtg }_ncoda_MERinl 
vvcur i_$ { idtg }_ncoda_MERinl 
laypr s_$ { idtg }_ncoda_MERinl 
subregional . depth .MERinl . a 
seaice_$ { idtg }_ncoda_MERinl 
zi42 . txt 

3949  'ilstn  '  =  i-origin  of  ncoda  subregion 

360  'jlstn  '  =  j-origin  of  ncoda  subregion 

1313  'idmn  '  =  i-extent  of  ncoda  subregion  (<=idm;  0  implies  idm) 

1569  'jdmn  '  =  j-extent  of  ncoda  subregion  (<=jdm;  0  implies  jdm) 

42  'kncoda'  =  number  of  ncoda  levels 
seatmp_$ { idtg }_ncoda_MERatl 
salint_$ { idtg }_ncoda_MERatl 
uucur i_${ idtg}_ncoda_MERatl 
vvcur i_$ { idtg }_ncoda_MERatl 
laypr s_${ idtg }_ncoda_MERatl 
subregional . depth .MERatl . a 
seaice_${ idtg}_ncoda_MERatl 
zi42 . txt 

2349  'ilstn  '  =  i-origin  of  ncoda  subregion 

360  'jlstn  '  =  j-origin  of  ncoda  subregion 

1751  'idmn  '  =  i-extent  of  ncoda  subregion  (<=idm;  0  implies  idm) 

1841  'jdmn  '  =  j-extent  of  ncoda  subregion  (<=jdm;  0  implies  jdm) 

42  'kncoda'  =  number  of  ncoda  levels 
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seatmp_$ { idtg }_ncoda_MERpal 
salint_$ { idtg }_ncoda_MERpal 
uucuri_${ idtg}_ncoda_MERpal 
vvcur i_$ { idtg }_ncoda_MERpal 
layprs_${ idtg }_ncoda_MERpal 
subregional . depth .MERpal . a 
seaice_${ idtg}_ncoda_MERpal 
zi42 . txt 

199  ' ilstn  '  =  i-origin  of  ncoda  subregion 
360  ' jlstn  '  =  j -origin  of  ncoda  subregion 

2525  ' idmn  '  =  i-extent  of  ncoda  subregion  (<=idm;  0  implies  idm) 

1841  'jdmn  '  =  j-extent  of  ncoda  subregion  (<=jdm;  0  implies  jdm) 

42  ' kncoda '  =  number  of  ncoda  levels 

seatmp_${ idtg}_ncoda_ARCatl 
salint_${ idtg}_ncoda_ARCatl 
uucuri_${ idtg}_ncoda_ARCatl 
vvcur i_$ { idtg }_ncoda_ARCatl 
layprs_${ idtg }_ncoda_ARCatl 
subregional . depth . ARCatl . a 
seaice_${ idtg}_ncoda_ARCatl 
zi42 . txt 


2440 

'  ilstn  ' 

=  i-origin 

of 

ncoda 

subregion 

2150 

'jlstn  ' 

=  j-origin 

of 

ncoda 

subregion 

1490 

'  idmn  ' 

=  i-extent 

of 

ncoda 

subregion 

(<=idm; 

0 

implies 

idm) 

551 

'  j  dmn  ' 

=  j-extent 

of 

ncoda 

subregion 

(<=j  dm; 

0 

implies 

jdm) 

42 

'  kncoda ' 

=  number 

of 

ncoda 

levels 

seatmp_${ idtg}_ncoda_ARCpac 
salint_${ idtg}_ncoda_ARCpac 
uucuri_${ idtg}_ncoda_ARCpac 
vvcur i_${ idtg}_ncoda_ARCpac 
layprs_${ idtg }_ncoda_ARCpac 
subregional . depth . ARCpac . a 
seaice_${ idtg}_ncoda_ARCpac 
zi42 . txt 

730  'ilstn  '  =  i-origin  of  ncoda  subregion 


2150  'jlstn  ' 

=  j-origin 

of 

ncoda 

subregion 

1335  'idmn  ' 

=  i-extent 

of 

ncoda 

subregion 

(<=idm; 

0 

implies 

idm) 

551  'jdmn  ' 

=  j-extent 

of 

ncoda 

subregion 

(<=j  dm; 

0 

implies 

jdm) 

42  'kncoda' 

=  number 

of 

ncoda 

levels 

seatmp  ${idtg} 

ncoda  ARCocn 

salint  ${idtg} 

ncoda  ARCocn 

uucuri  ${idtg} 

ncoda  ARCocn 

vvcur i  ${idtg} 

ncoda  ARCocn 

layprs  ${idtg} 

ncoda  ARCocn 

subregional . depth . ARCocn . 

a 

seaice  ${idtg} 

ncoda  ARCocn 

zi42 . txt 

2459  'ilstn  ' 

=  i-origin 

of 

ncoda 

subregion 

2650  'jlstn  ' 

=  j-origin 

of 

ncoda 

subregion 

1600  'idmn  ' 

=  i-extent 

of 

ncoda 

subregion 

(<=idm; 

0 

implies 

idm) 

1297  'jdmn  ' 

=  j-extent 

of 

ncoda 

subregion 

(<=j  dm; 

0 

implies 

jdm) 

42  'kncoda' 

=  number 

of 

ncoda 

levels 

NONE 

E-o-D 

/bin/mv  $ {arc2 } . $ { idtg} .a  $ {arc2 } . $ { idtgtodpl } $ {HR}_$ { idtg} $ { HR} .a 

/bin/mv  $ {arc2 } . $ { idtg} .b  $ {arc2 } . $ { idtgtodpl } $ { HR } _ $ { idtg} $ { HR} . b 

# 

date 

# 

#  make  links  from  files  like  archv . 2002 01 01 . a  to  archv. 2002_001_00 . a 

# 

cd  $ { INP } 

# 


64 


NRL/MR/7320— 09-9175 


GOFS  V3.0  User’s  Manual 


foreach  FILE  ('Is  $ {arc2 } . $ { idtgtodpl } $ {HR}_$ { idtg} * ' ) 
set  iyear='echo  ${FILE}  | awk  -F${HR}_  1 { a=substr ( $NF, 1 , 4 ) ; print  a}1' 
set  ii='echo  ${FILE}  | awk  -F${HR}_  1 { a=substr ( $NF, 5, 2 ) ; print  a}1' 
set  imon='echo  $ { i i }  |  sed  1 s/A0*\ ( . *\) $/\l/ ' ' 

set  ii='echo  ${FILE}  | awk  -F${HR}_  1 { a=substr ( $NF, 7 , 2 ) ; print  a}1' 

set  iday='echo  $ { i i }  |  sed  1 s/A0*\ ( . *\)  $/\l/ '  ' 

set  idtg= ' $ { BIN } / ymd2doy . csh  ${iyear}  ${imon}  ${iday}' 

set  ab='echo  ${FILE}  |  sed  -n  ' s/A . *\ . \ ( . *\) $/\l/p ' ' 

echo  ${iyear}  ${imon}  ${iday}  ${ab}  ${idtg} 

/bin/rm  . . /data/ incup/ incupd. $ { idtg}_$ {HR} . $ {ab} 

In  -s  ../../${ case }/${ FILE }  . . /data/incup/incupd. $ { idtg}_$ { HR} . $ { ab} 
end 
# 

date 

# 

# 

#  don't  use  this  part  unless  you  have  dedicated  time  and  you  want  to 

#  maximize  the  the  use  of  the  remainder  of  the  time 
#if(${idtgl}  ==  ${ idtgtodpl }) then 

#echo  '****************  FORECAST  **********************' 

#echo  ' idtg=  '${idtg}  '  >  idtgtod=  '  ${idtgtod} 

## 

#set  targhr=ll 
#set  targmin=60 
#set  nowhr='date  +%H' 

#set  diffhr='expr  ${targhr}  -  ${nowhr}' 

#set  diffhr='expr  ${diffhr}  -  1' 

#set  nowmin='date  +%M' 

#set  dif fmin= ' expr  ${targmin}  -  ${nowmin}' 

#set  dif fmin= ' expr  ${diffmin}  -  5' 

#if (  ${diffmin}  <  0  )  then 

#  set  dif fmin= ' expr  60  -  ${diffmin}' 

#  set  diffhr='expr  ${diffhr}  -  1' 

#endif 

## 

#set  dif fmin= ' echo  ${diffmin}  |  awk  1 {printf  "%02d  \n",$l}’' 

#set  runtim=$ { dif fhr } : $ { dif fmin } 

#echo  ${runtim} 

##number  of  forecast  days  of  model  integration 
#setenv  frcstdays  2 
##setenv  frcstdays  3 
## 

##if ( $ { dif fhr }  >=  4  &&  ${diffmin}  >=  45)  then 

##  setenv  frcstdays  4 

##endif 

##if ( $ { dif fhr }  ==  5  &&  ${diffmin}  >=  00)  then 

##  setenv  frcstdays  4 

##endif 

##if ( $ { dif fhr }  >=  6  &&  ${diffmin}  >=  20)  then 

##  setenv  frcstdays  5 

##endif 

#date 

#echo  diffhr  ${diffhr}  diffmin  ${diffmin}  runtim  ${runtim}  frcstdays  ${frcstdays} 

#  /bin/rm  $ {LOGS } /$ { E }pbs_ncoda_$ { idtgtodpl } $ {HR}_$ { idtgl } $ {HR} * . { com, log} 

##  awk  -f  $ {MVOI }/.. /HYCOM. awk_new  runtim=$ { runtim}  hr=${HR}  nmdays=$ { frcstdays } 
t0=${E}pbs_${idtgl}  tod=$ { idtgtodpl }  tl=${ idtgl}  $ {MVOI }/../$ {E }pbs_ncoda . com  > 

$ {LOGS } / $ { E } pbs_ncoda_$ { idtgtodpl } $ {HR}_$ { idtgl } $ {HR} . com 

#  awk  -f  $ {MVOI }/.. /HYCOM . awk  runtim=$ { runtim}  hr=${HR}  nmdays=$ { frcstdays } 
t0=${E}pbs_${idtgl}  tod=$ { idtgtodpl }  tl=${ idtgl}  $ {MVOI }/../$ {E }pbs_ncoda . com  > 

$ {LOGS } /$ { E } pbs_ncoda_$ { idtgtodpl } $ {HR}_$ { idtgl } $ {HR} . com 

##FORECAST 

#endif 

# 
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#  submit  hycom  for  next  day 

# 

cd  $ { LOGS } 

#OMS 

/u/home/wallcraf /bin/q_navo  $ { E }pbs_ncoda_$ { idtgtodpl } $ {HR}_$ { idtgl } $ {HR} . com 

# 

#  submit  the  plotting  job 

#  hindcast 

/u/home/wallcraf /bin/ q_navo  $ { E }plotpost_$ { idtgtodplt } $ { HR } _ $ { idtgl } 00 . csh 

#  real  time 

# /u/home/wallcraf /bin /q_navo  $ { E } plotpost_$ { idtgtodpl } $ { HR } _ $ { idtgl } 00 . csh 
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APPENDIX  F 

Script  727pbs_ncoda_  ${idtgtod}  1 8  $ {idtg}  18.com 

This  is  a  Script  to  run  HYCOM. 

#!  /bin/csh  -x 
#PBS  -N  999pbs 
#PBS  -j  oe 

#PBS  -o  999pbsXX.log 
#PBS  -1  mppwidth=619 
#PBS  -1  mppnppn=8 
#PBS  -1  walltime=2 : 00 : 00 
#PBS  -W  umask=027 
#PBS  -A  NRLSS018 
#PBS  -q  internal3d 
# 

#BSUB  -J  727ncoda 
#BSUB  -M  1500000 
#BSUB  -a  poe 
#BSUB  -n  379 

#BSUB  -R  "rusage [ntbl_windows=32 ] span [ptile=l 6 ] " 

#BSUB  -W  4:00 

#BSUB  -q  internal3d 

#BSUB  -P  NAVOSHYC 

# 

set  echo 
set  timestamp 
date 
C 

C  -  Preamble. 

C 

setenv  OS  'uname' 
switch  ($OS) 
case  1 AIX 1 : 
hostname 

#  setenv  TMPDIR  /scr/${user} 
breaksw 

case  1  Linux ' : 
which  yod 

if  (!  $status)  then 
setenv  OS  XT3 

#  setenv  TMPDIR  /tmp 

else 

#  setenv  TMPDIR  /tmp 

endif 

which  aprun 

if  (!  $status)  then 

#  setenv  OS  XT4 
setenv  OS  XT5 

#  setenv  TMPDIR  /scr 

endif 

breaksw 
default : 

echo  'Unknown  Operating  System:  '  $OS 
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echo  'configured  for  AIX  only' 
exit  (1) 
endsw 


C 

if  ($?J0BNAME)  then 

setenv  PBS_JOBNAME  $ { JOBNAME } 
setenv  PBS_JOBID  $$ 
endif 

echo  PBS_ JOBNAME  $PBS_JOBNAME  PBS_JOBID  $PBS_JOBID 
C 

C  -  Automatic  Run  Script. 

C  -  Submit  with  msub,  or  msub_csh,  or  msub_ll  command. 

C  -  Multiple  segment  version,  set  no.  of  segments  on  foreach  below. 

C 

C  -  Set  up  for  incremental  updating 

C 

C 

C  -  E  is  expt,  P  is  permenant  directory,  S  is  /trap  directory. 

C 

# 

setenv  E  727 
setenv  REG  GLBa0.08 
setenv  EXPT  expt_72.7 
# 

#  -  set  by  awk  script  for  daily  run 

# 

setenv  nmdays  1 
setenv  idtg  20031102 
setenv  idtgtod  20031102 
setenv  HR  18 
setenv  HR2  00 


# 

#  -  restart  is  "tomorrow"  since  we  start  at  18Z 

# 

set  idtgrestart=' /u/home/$ { user } /bin/addndays  yyyymmdd  ${idtg}  tl ' 
echo  idtg  ${idtg}  idtgtod  ${idtgtod}  idtgrestart  $ { idtgrestart } 

# 

set  LOGS=/scr/$ {user} /hycom/$ {REG} /$ { EXPT } /logs 
mkdir  -p  ${LOGS} 

set  RUN=/u/home/$ {user} /hycom/ $ {REG} / $ { EXPT } 

# 

#  -  remove  temporary  NCODA  files 

# 

Cd  $ { RUN } 

/u/home/wallcraf /bin/q_navo  $ {E } removef iles_2 . com 

# 


setenv  P  $cwd 
switch  ($OS) 
case  'AIX': 
case  ' XT5 ' : 

# 

setenv  S  'echo  $cwd  | 
breaksw 
e  1  XT  3  ’  : 
case  ' XT4  '  : 
case  ' OSF1 ' : 

# 

setenv  S  'echo  $cwd  | 


substitute  /scr  for  /u/home 
awk  '{print  "/scr"  substr ($0, 8, length) }' ' 


substitute  /work  for  /??????? 
awk  '{print  "/work"  substr ( $0 , 3 , length)  }' 
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breaksw 
case  ' IRIX64 ' : 

#  substitute  /scr  for  /u/home 

setenv  S  'echo  $cwd  |  awk  '{print  "/scr"  substr ($0, 8, length) }' ' 
breaksw 
case  1 unicos 1 : 

#  substitute  /trap  for  /u/b 

setenv  S  'echo  $cwd  |  awk  '{print  "/trap"  substr ($0, 5, length) }' ' 
breaksw 
endsw 
C 

Is  -laFq 
C 

C  -  check  the  RUNNING  flag. 

C 

if  (  -e  RUNNING  &&  !  -e  RUNNING_$PBS_JOBID  )  then 
C 

C  -  MODEL  IS  ALREADY  RUNNING  -  EXIT. 

C 

exit 

endif 

touch  RUNNING 

touch  RUNNING_$PBS_JOBID 

C 

C  -  Generate  the  next  model  script. 

C 

setenv  Y01  103 
setenv  AB  a 

# 

setenv  SCRIPT  $ { E } y$ { idtg } . com 
/bin/rm  -f  ${SCRIPT} 

awk  -f  $ { RUN} /$ {E } . awk  nmdays=$ { nmdays }  hr=${HR}  y01=${Y01}  tod=$ { idtgtod} 
ab=${AB}  td=$ { idtg}  ${RUN}/${E} .com  >  ${SCRIPT} 

# 

#  -  Run  the  Script. 

# 

set  script  =  $SCRIPT 

set  reqname  =  $ { PBS_JOBNAME } 

In  -fs  ${ reqname }. log  $script : r . log 

# 

mkdir  -p  $S 

cp  $ { SCRIPT }  $S/ $ { SCRIPT } 
cd  $S 
# 

# - 

csh  $ { SCRIPT } 

# - 

# 

wait 

date 

# 

#  -  submit  next  day  of  todays  run 

# 

#  -  submit  next  day  of  todays  run 

#  -  update  only  idtg  if  a  real  time  run  that  goes  back  #  of  days 

# 

set  idtg= ' /u/home/$ {user } /bin/addndays  yyyymmdd  ${idtg}  +1' 

# 
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echo  '**************  in  lsf_ncoda.com  ****************' 

# 

echo  idtg  ${idtg}  shellsub_hindcast . csh 
$ {RUN } /mvoi/ $ {E } shellsub_hindcast .csh  $ { idtg} 

# 

#echo  idtg  ${idtg}  idtgtod  ${idtgtod}  shellsub_daily . csh 
#$ { RUN } /mvoi/$ { E } shell sub_daily .csh  ${idtg}  ${idtgtod} 

# 

cd  $P 

# 

Is  -laFq 

# 

#  Final  Clean  Up. 

# 

/bin/rm  -f  RUNNING 

/bin/rm  -f  RUNNING_$PBS_JOBID 

# 

#  Exit. 

# 

exit 
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